arXiv:1501.07266vl [cond-mat.soft] 28 Jan 2015 


A corresponding version of this manuscript has been accepted for publication by Physics Reports. 
Journal Reference: Phys. Rep. 554 , 1 (2015) 

URL: http://www.sciencedirect.com/science/article/pii/S0370157314003871 
DOI: 10.1016/j.physrep. 2014-10.001 


Tuned, driven, and active soft matter 


Andreas M. Menzel 

Institut fur Theoretische Physik II: Weiche Materie, Heinrich-Heine-Universitdt Diisseldorf, 
Universitatsstrafie 1, D-40225 Diisseldorf, Germany; 
phone number: +49-211-81-12056 


Abstract 

One characteristic feature of soft matter systems is their strong response to external stimuli. 
As a consequence they are comparatively easily driven out of their ground state and out of 
equilibrium, which leads to many of their fascinating properties. Here, we review illustrative 
examples. This review is structured by an increasing distance from the equilibrium ground 
state. On each level, examples of increasing degree of complexity are considered. In detail, 
we first consider systems that are quasi-statically tuned or switched to a new state by 
applying external fields. These are common liquid crystals, liquid crystalline elastomers, 
or ferrogels and magnetic elastomers. Next, we concentrate on systems steadily driven 
from outside e.g. by an imposed flow field. In our case, we review the reaction of nematic 
liquid crystals, of bulk-filling periodically modulated structures such as block copolymers, 
and of localized vesicular objects to an imposed shear flow. Finally, we focus on systems 
that are “active” and “self-driven” . Here our range spans from idealized self-propelled 
point particles, via sterically interacting particles like granular hoppers, via microswimmers 
such as self-phoretically driven artificial Janus particles or biological microorganisms, via 
deformable self-propelled particles like droplets, up to the collective behavior of insects, fish, 
and birds. As we emphasize, similarities emerge in the features and behavior of systems 
that at first glance may not necessarily appear related. We thus hope that our overview will 
further stimulate the search for basic unifying principles underlying the physics of these soft 
materials out of their equilibrium ground state. 
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1. Introduction 


It is easy to motivate the study of soft matter systems. Everyday life provides us with 
the most immediate examples. When we go shopping, we carry home our groceries in plastic 
bags made from polymers. Most flat screen displays exploit the switching behavior of liquid 
crystals in an electric field. We use solutions of surfactants when we soap ourselves while 
taking a shower. And finally, a major part of ourselves is actually soft matter, as almost all 
biological cells are. Considering the shear amount of soft matter distributed commercially, 
research activities in this field naturally extend far beyond purely academic purposes. 

It is much less straightforward to give a definite, unambiguous, and precise definition of 
what the term “soft matter” actually means or which states of matter it actually includes. 
Maybe one should proceed the other way around and explain which materials it excludes 
by comparison to common other states of matter. 

On the one hand, we interpret the term “softness” as a strong response to a comparatively 
weak force |lj. This becomes most illustrative when we compare the elastic behavior of metal 
solids and ordinary rubbery materials. The Young modulus of a rubber can for example 
be ten thousand times lower than that of steel |2j. I.e. applying a tensile force to that 
piece of rubber can stretch it about ten thousand times further than when applying it to a 
corresponding piece of steel. In general, due to these huge responses to external forces, the 
observed behavior of soft materials is markedly nonlinear [I]. 

On the other hand, the response of soft matter systems to external stimuli is generally 
comparatively slow pQ. The relaxation times in polymers can be hours or days (or even 
more) [[3] , whereas it takes nanoseconds in simple liquids pfj. It is therefore easy to drive 
and maintain the systems out of equilibrium, and non-equilibrium effects play a central role 
in the study of soft materials (Tj. 

One reason for these properties is the size of the building blocks [TJ. These sizes typically 
fall into the range of several nanometers up to about a hundred micrometers. Some examples 
together with a rough classification of the possible associated dimensions are schematically 
depicted in Fig. [lj Clearly, small rod-like liquid crystalline molecules like PAA with a 
length of about two nanometers are on the lower end [5]. Colloids, i.e. mesoscopic particles 
stabilized in a permanent dispersion, cover a broad range of the considered sizes by definition 
M- On the contrary, a rubber made of sufficiently chemically crosslinked polymer chains 
can actually be considered as one giant molecule and reach macroscopic dimensions [3]. At 
this point a strict definition via the size of the constituents becomes more involved. 

Nevertheless, we can often say that soft matter building blocks are large enough so 
that quantum effects can be neglected, but small enough that thermal fluctuations play 
a significant role pm isg. For example, most of the elastic behavior of rubbers must be 
attributed to entropy, i.e. thermal fluctuations M- Entropy also supports the stabilization 
of colloidal dispersions as a complement to the necessary artificial stabilization mechanisms 

M- 

Apart from that, the size of the building blocks allows us to use coarse-grained descrip¬ 
tions to characterize their collective physical behavior in. Chemical details can often be 
neglected, which leads to a certain degree of “universality” in the characterization. This 
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Figure 1: Schematic illustration of typical examples of soft-matter building blocks and materials, together 
with a rough estimate of their characteristic sizes (horizontal bars indicate a broader range of possible 
occurring length scales). The rod-like molecule PAA is a common representative of low-molecular-weight 
liquid crystals 0 that can show an orientationally ordered anisotropic nematic phase. Surfactant molecules 
combine two antagonistic features on one building block: their head is usually hydrophilic, whereas their 
tail(s) is (are) mostly hydrophobic [5]- This leads, for instance, to the formation of closed bilayer membranes 
in the form of vesicles uni, with giant unilamellar vesicles (GUVs) as an extreme example (in reality, unlike 
the schematic, the layer thickness is orders of magnitude smaller than the radius) [TTtiT3] . Polymers [JMl] in 
the simplest case can be treated as flexible linear chain-like objects. When covalently chemically crosslinked 
in the form of rubbers and elastomers (crosslinks are indicated by dots) a macronrolecule of macroscopic 
dimension can be obtained M- Block copolymers feature chemically different parts (marked by different 
colors) that can micro-phase separate into regularly ordered structures |15| , here indicated for a lamellar 
texture. In colloids, the constituting particles are stabilized to form a permanent dispersion IMj. Again, 
different sides of the particles can be functionalized with different chemical properties so that the resulting 
Janus particles [16] form clusters or even regular lattices mm- Finally, biological cells such as motile 
bacteria EH2U combine a multitude of functionalities. Collective arrangements of biological cells in the 
form of biological tissue reaches macroscopic dimensions. 
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point of view got particularly famous in the field of polymer physics [23]. 

Finally, the size of the constituents makes it possible to arrange different or even an¬ 
tagonistic properties on one building block. Some examples are depicted in Fig. [l] These 
are block copolymers, where chemically different connected parts of the molecules tend to 
micro-phase separate into regular spatial textures [13]; surfactant molecules composed of 
hydrophobic (“water-fearing”) and hydrophilic (“water-loving”) parts that self-assemble for 
instance into bilayer membranes and vesicles in an aqueous environment PUS]; and colloidal 
Janus particles [16] that can form finite clusters or regular lattices mm- 

With all these facts in mind, it is not surprising that the rheological behavior of such 
systems is relatively complex when compared to simple liquids [21]. Many soft matter 
systems are therefore also classified as “complex fluids”. Sometimes the terms “soft matter” 
and “complex fluids” are even used interchangeably [23]. 

In the following, we will concentrate on soft matter systems out of their equilibrium 
ground state. It is impossible to cover all aspects and areas of this subfield, which was 
expanding so rapidly over the past decades. Instead, we will mainly focus on typical rep¬ 
resentatives in each part and embed these examples into the broader context. Roughly, 
we will increase the degree of activation that leads to non-equilibrium states through the 
subsequent sections. 

First, we consider the switching or tuning of the states of soft matter systems by external 
fields. For example, as noted above, in display devices the orientation of the optical axis 
in liquid crystals is switched by external electric fields. Through a static external held, the 
system is forced out of its initial equilibrium ground state. However, it can reach a new 
static equilibrium state in the presence of the external held. We will mainly concentrate 
on liquid crystalline as well as superparamagnetic or ferromagnetic gels and elastomers, the 
optical and mechanical properties of which can be tuned by external helds. They can serve 
as soft actuators, and even their possible role as “artificial muscles” was pointed out. 

Next, we focus on systems permanently driven out of equilibrium by steady external shear 
hows. A brief account of the rheological behavior of liquid crystals and block copolymers 
is included. Coming back to the introductory cases mentioned above, an example situation 
from everyday life would be soap or body lotion while we are spreading it on our skin. 
Related to this topic, we address the progress made on the dynamics of vesicles when driven 
by an external shear how. 

The broadest section is devoted to the recently exploding held of “active” soft matter 
systems. We understand the word “active” in the sense that individual building blocks have 
a sort of propulsion mechanism, which leads to motion. Crawling cells or swimming bacteria 
can run this mechanism for instance by consuming fuel or food, whereas colloidal particles 
can self-propel via self-phoretic effects. In general, a key point is that the direction of motion 
is not prescribed from outside. These ingredients lead to new and fascinating properties of 
the individual constituents and when they collectively act together. There is a significant 
overlap between the helds of soft matter and biological physics in this area. 

Finally, a short summary is appended. We conclude by giving a brief impression of the 
current state of the held and its evolution. 
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2. Timed soft matter 


By “tuning” we understand an adjustment of the material properties to a state that 
serves the current purpose or convenience. The crucial point with the soft matter systems 
introduced below comes with the non-invasive way of this adjustment. An external held 
is applied to the closed cell or piece of material to achieve the new setting. In general, a 
static external held can force the system to acquire a new static equilibrium state in the 
presence of the held. If the external held is increased or decreased in small steps, after each 
adjustment the system hnds itself in a weakly non-equilibrium situation. It has to find the 
new equilibrium state in the presence of the switched external held. Frequently, this process 
can be described by simple relaxation dynamics following the minimum of the modified free 
energy landscape. If the external held is adjusted slowly enough, a quasi-static switching 
can often be achieved, avoiding genuinely non-equilibrium situations. If our focus is on the 
process of conversion of external held energy to reach a new state of the system, we may 
speak of “actuation”. 

2.1. Low-molecular-weight liquid crystals 

As a starting point, let us briehy recall the well-known Freedericksz transition in nematic 
liquid crystals [5j. This transition can be induced from outside, for instance by applying 
an external electric held. As will be explained in more detail below, it involves a collective 
reorientation on the molecular level, changing the optical birefringence properties of a sample 
cell. In this way, the light transmittance of the cell can be reversibly tuned from outside. 

In the most illustrative case, the liquid crystalline building blocks are small rod-like 
organic molecules. The molecular structure of one characteristic example, namely PAA, 
was depicted in Fig. [l] It shows the typical dimensions, which in this case are about two 
nanometers in length and about half a nanometer in thickness p3j- At high enough temper¬ 
atures, the material is in an isotropic disordered liquid state as schematically illustrated in 
Fig. @ (a). The system does not feature any long-ranged positional nor orientational order of 
the constituents. It behaves in an isotropic way. In contrast to that, we find both positional 
and orientational order at low enough temperatures when the material forms a crystal. Then 
the molecules are stacked on regular lattice sites and orient their long axis along a common 
direction, see Fig. [2] (d). Between these two extremes of a disordered liquid and a crystalline 
solid, liquid crystals feature further intermediate phases, the so-called mesophases. These 
liquid crystalline phases show a higher degree of order (corresponding to a lower degree of 
symmetry) than the fluid but a lower degree of order (corresponding to a higher degree 
of symmetry) than the crystal. For example, in the nematic state, the one that we will 
refer to in the following, there is no long-ranged positional order. However, on average, the 
molecules orient along a common direction, see Fig. [2] (b). In this way, the material becomes 
anisotropic. Smectic states, which we will come back to later in this review, feature a regular 
stacking in layers, see Fig. [2] (c). More precisely, they show a quasi-long-ranged positional 
order in one spatial dimension. 

When the system is in the nematic phase, the building blocks tend to collectively order 
their long axes along a common direction. This average orientation is called the nematic 
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Figure 2: Schematic illustration of typical phases of common low-molecular-weight liquid crystals consisting 
of elongated rod-like molecules, (a) The disordered isotropic liquid-like state shows neither long-ranged 
orientational nor long-ranged positional order, (b) In the nematic phase, there is no long-ranged positional 
order. However, apart from thermal fluctuations, the long axes of the molecules on average orient along 
a common direction called the director, (c) Smectic phases are characterized by a layer-like positional 
arrangement of the molecules, while in the depicted smectic-C state there is no long-ranged positional order 
within each layer, (d) Finally, crystalline solids show the highest degree of order with long-ranged positional 
order in the form of stacking of the molecules on regular lattice sites and long-ranged orientational order. 
Transitions between the depicted phases can be induced in corresponding liquid crystalline materials for 
example by temperature changes. 


director n [5|. Genuinely nematic phases are non-polar, that is the directions +n and —n 
cannot be distinguished. The degree of orientational order is measured by a scalar order 
parameter 

s = ^((3cos 2 tf-l)), (1) 

where the average (...) is taken over all molecules in the sample, and d for each individual 
molecule measures the angle between its long axis and the average orientation n. We obtain 
s = 0 in the disordered isotropic liquid state, and s = 1 in a perfectly ordered nematic 
state. A combined symmetric traceless order parameter tensor S can be introduced, which 
contains both information, the degree of orientational order and the average orientation n, 



Here, nn denotes a dyadic product, d is the dimension of space (usually d = 3), and I 
represents the unity matrix^] 

On the one hand, since n corresponds to the average orientation of the ordered molecules, 
it is an axis of anisotropy. For example, the dielectric constants eii and e± measured parallel 
and perpendicular to n, respectively, generally differ from each other. Let us consider the 
case in which the dielectric anisotropy e a = en — e_L is positive, e a > 0. Then a static external 


^lso another notation is common, in which the letters S and Q are used instead of s and S, respectively. 
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Figure 3: Reorientation of the nematic director induced by an external electric field in a thin cell of prescribed 
boundary conditions. The nematic director is indicated by the double-headed arrow. At the upper and lower 
confining plate it is rigidly anchored horizontally. On the one hand, if the sample is thin enough, this leads 
to a homogeneous horizontal alignment of the director to avoid spatial distortions in the director orientations 
(left). On the other hand, for positive dielectric anisotropy of the nematic molecules, a vertically applied 
external electric field E tends to orient the nematic molecules and the director in vertical direction (center). 
Both effects - horizontal alignment due to the boundary conditions and vertical alignment through the 
external electric field - compete with each other. Balancing the two resulting torques leads to a finite 
angle of director reorientation with respect to the horizontal direction (right) above a non-vanishing critical 
threshold amplitude of the electric field. 


electric field E tends to orient the director n along itself. This tendency is reflected by a 
corresponding contribution to the effective thermodynamic potential 

Fei e c = ~^e a (n ■ E) 2 . (3) 

In this expression, possible constant prefactors due to the system of measure are absorbed 
into the dielectric anisotropy e a . 

On the other hand, spatial variations of the director, i.e. spatial distortions, increase 
the thermodynamic potential and are therefore suppressed. One way of quantifying this 
thermodynamic penalty is by contracting the tensor Vn, 

Fdist = \K{V n):(Vn), (4) 

with K > 0 another material-dependent parameteij^J As an example, the nematic liquid 
crystal can be confined in a thin cell between two parallel plates as depicted in Fig. [3j The 
surfaces of the cell can be prepared such that they anchor the orientation of the director, 
for instance in a parallel way, as indicated on the left-hand side of Fig. [3} This, via Eq. Q, 
sets a uniform director orientation in the whole cell. 

If now an external electric field is switched on with an orientation perpendicular to 
the cell surfaces, see the center of Fig. [3j it tends to reorient the director and align it 
parallel to itself as given by Eq. (|3]) . However, via Eq. Q, such a reorientation induces 


2 We here use a simplified version of this thermodynamic expression. Due to the anisotropy of the system, 
the coefficient K in general is replaced by a fourth-rank tensor [5]. 
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thermodynamically penalized distortions because of the surface anchoring on the plates. 
The competition between the two effects leads to a critical electric field amplitude that has 
to be exceeded for director reorientations to become possible. Above this threshold field 
amplitude, director reorientations are induced. From Eqs. (J3| and Q this critical field 
amplitude follows directly as 


E c = 


3 #- 

d\e a 


(5) 


where d is the thickness of the cell and the impact of possible electric currents has been 
neglected. In the new static reoriented state, the two effects of the reorienting electric held 
and the counteracting surface anchoring balance each other. This is indicated on the right- 
hand side of Fig. [3j When the electric held is switched oh, the director rotates back into its 
initial state due to the surface anchoring. 

In this way, the anisotropy properties of the presented sample cell can be reversibly 
tuned from outside. Naturally, this anisotropy in a real sample is connected to optical 
birefringence. Thus the optical transmission properties of the cell can be reversibly switched 
and adjusted by the external electric held. This is the underlying principle extensively 
used to construct optical display devices that are based on liquid crystals [26], However, 
the technical realizations today involve much more complicated geometries than the simple 
Freedericksz cell outlined above. 


2.2. Liquid crystalline elastomers 

Commonly used liquid crystals have relatively low molecular weight and are classihed as 
complex huids. As suggested by the name, they typically how away when not maintained by 
a containing cell. Another promising material would be one that is stable and self-standing 
without supporting container walls, but still features the optical transmission properties of 
common liquid crystals. Such materials were actually realized by stepping towards higher 
molecular weights in the form of liquid crystalline polymers and elastomers. 

In general, polymers are synthesized by covalently binding together from hundreds up 
to millions identical molecular repeat units, the so-called monomers. In this way, new huge 
macromolecules are formed [3]. Different chemical routes for this procedure are available. 
The number of monomers on the resulting polymeric molecule defines the degree of poly¬ 
merization. In the simplest case, long linear chain-like molecules are obtained, but several 
further molecular architectures can be realized. 

The simplest theoretical model describing linear polymer chains is the freely-jointed chain 
model [3], It already gives us enough qualitative insight to understand the discussion below. 
In this model, as illustrated in Fig. [4j the polymeric molecule is reduced to a chain of N 
straight segments that arc linked to each other at their ends. Each segment has an identical 
length b much larger than the size of a single monomer. We denote the positions of the joints 
between the segments as R„ for n = 1,..., N — 1, while R 0 and Rat give the positions of the 
two ends of the whole resulting chain, see Fig. [4j The chain is “freely” jointed in the sense 
that each segment orientation is independent of the orientation of all other segments. Thus, 
for each segment, the probability distribution i/j(r n ) of its end-to-end vector r n = R n — R n _i 
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Figure 4: Illustration of the freely-jointed chain model for a linear polymer chain. The chain is represented 
by N freely-jointed segments of end-to-end vectors r„ = R„ — R„_i (n = 1 where R„ (n = 

1 ,N — 1) give the positions of the joints. All segments are of identical length, ||r„|| = b for all n, and 
their orientations are completely independent of each other. R = R^v — Ro is the end-to-end vector of the 
whole polymer chain. 


is given by 

i>( r n) = ^pH\\r n \\-b) (n = 1 ,..., N), ( 6 ) 

where £(■) denotes the Dirac ^-function. It serves here to prescribe the length b of the 
segment. In this minimum model, the probability distribution to find a certain chain con¬ 
figuration is simply given by the product of these individual single-segment probability 
distributions i/j( r n ). 

On this basis, it is straightforward to determine the probability distribution for the 
end-to-end vector R = — Rq of the whole chain. For IV 1, we obtain 


^(R) 


/ 3 \ 3/2 

(j 2 exp 


3R 2 1 

2Nb 2 J ' 


(7) 


Obviously, this is of Gaussian form. 

From a statistical mechanics point of view, we expect the probability distribution to 
be of the form -0(R) ~ exp[— F(TV)/k B T\. Here F( R) denotes the energy of the chain 
configuration with end-to-end vector R, ks is Boltzmann’s constant, and T the temperature. 
Comparison with Eq. (J7]) shows that 


F(R) 


1 3k B T 2 
2 Nb 2 


( 8 ) 


Thus the polymer chain behaves like an elastic spring between its two end points Ro and 
Ryv- In the unstrained state the length of this spring vanishes. We find 3k B T/Nb 2 for the 
magnitude of its elastic constant. 

Based on this simple model, we see that polymeric materials show an clastic response. 
The important point is that these elastic forces mainly do not origin from an internal energy. 
They are mostly entropic in nature. Everything that enters Eq. (|7|) is the entropy resulting 
from the single-chain configuration. The elastic constant in Eq. (ISJ) is set by the temperature 
T. 
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Figure 5: Schematic illustration of two architectures of liquid crystalline polymers. Again the liquid crys¬ 
talline rod-like components are indicated by ellipsoids as in Fig. [2] Non-liquid-crystalline parts of the 
polymer “backbone” chains are depicted by the bright strings, (a) On the one hand, in main-chain liquid 
crystalline polymers, the liquid crystalline units are part of the polymer backbones, (b) On the other hand, 
the liquid crystalline units are attached as side groups to the polymer chains in side-chain liquid crystalline 
polymers. Usually this is achieved via shorter hydrocarbon chains, so-called spacer groups, which here are 
indicated in darker color. 


As a next step, one would like to obtain an elastic material of macroscopic dimension. 
For this purpose, many polymer chains are chemically crosslinked, i.e. covalently bound, to 
form one chemical unit. In this way, an elastomer or rubber is obtained. The main source 
of elasticity is again of entropic origin and referred to as “rubber elasticity” [2]. 

To make the connection to the held of liquid crystals, first liquid crystalline polymers 
were synthesized. Two principal architectures are available: on the one hand, liquid crys¬ 
talline molecules can be part of the polymer backbone chains in liquid crystalline main-chain 
polymers [27], see Fig. [5] (a); on the other hand, the liquid crystalline molecules can be at¬ 
tached as side-groups to a polymer backbone in liquid crystalline side-chain polymers j2S], see 
Fig. [5] (b). When in a later step such liquid crystalline polymers were chemically crosslinked, 
rubbery liquid crystalline elastomers were obtained [29, 30]. In this way, the properties of 
liquid crystals and those of elastic rubbers could be combined in one material. The coupling 
provides new interesting features as described below. 

There are different protocols available to obtain different kinds of liquid crystalline elas¬ 
tomers concerning the nematic orientational alignment in the samples. Without special 
aligning mechanisms it is natural for nematic bulk materials in general not to feature a 
globally aligned nematic director. Already for low-molecular-weight liquid crystals, bulk 
samples in the nematic phase usually appear turbid. They strongly scatter light [5]. The 
reason are thermal fluctuations in the director orientations. As a consequence, the nematic 
director is not macroscopically aligned, but continuously varies over space. Additional effort, 
such as the prescribed alignment on the sample surfaces depicted in Fig. [3] is necessary to 
achieve a global alignment of the director orientation. 

Naturally, thermal orientational fluctuations also occur in polymeric substances. There¬ 
fore, if during the crosslinking process of synthesizing a liquid crystalline elastomer no special 
action is taken, polydomain samples are obtained [3T|. The director in these polydomain 
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materials is not globally uniform in the nematic state but continuously varies over space. 
However, protocols to manufacture monodomain samples featuring a macroscopically aligned 
director are available. For this purpose, during the final crosslinking step of the manufactur¬ 
ing process, the sample must be maintained in the liquid crystalline phase and the director 
must be homogeneously aligned in the polymer solution. In practice, this was achieved 
by external magnetic fields \T2\ UFO . external electric fields [33], aligning boundary condi¬ 
tions [351I3T] . or stretching a pre-crosslinked sample [38] [39]. A corresponding memory of 
the director orientation during crosslinking is then imprinted into the material during the 
crosslinking procedure. When the liquid crystalline molecules were reoriented by some ex¬ 
ternal action, they were observed to rotate back into the memorized orientations after the 
external influence is switched off again. This was experimentally verified even for polydo¬ 
main samples crosslinked in the liquid crystalline state [HOIST]. Apparently the orientations 
are stored in the architecture of the polymer network. 

The most interesting properties of liquid crystalline elastomers arise from their orient- 
ational-deformational coupling. On the one hand, a mechanical deformation like stretching 
or compression can tune the state of birefringence and the optical properties. As one specific 
example, liquid crystalline elastomers were used to experimentally demonstrate the possibil¬ 
ity to construct mirrorless lasers of a wave-length that is tunable by mechanical deformations 
@2103]. More in general, it was shown for various different samples that mechanical de¬ 
formations can reorient the director [331 EZH32103H2S]. This is indicated on the left-hand 
side of Fig. [6] when a nematic liquid crystalline elastomer is stretched perpendicularly to its 
initial director orientation, a reorientation of the director into the stretching direction can 
result. The reorientation process involved a pronounced nonlinearity in the corresponding 
stress-strain curves I3ZM0Z]. On the other hand, it was demonstrated that the inverse ef¬ 
fect is also within reach experimentally: reorientations of the director can induce mechanical 
deformations. From an application point of view, this would allow the construction of soft 
actuators. One could use an external electric field to reorient the director, which is trans¬ 
formed into mechanical actuation. Unfortunately, so far the necessary field amplitudes were 
found to be too high for practical purposes. However, for samples swollen with a common 
low-molecular-weight (small-molecule) liquid crystal, a sizeable amount of so-called electro¬ 
mechanical coupling was observed [3~5l !48l 157] . In this case the director can be reoriented by 
an external electric held, which can result in deformations of up to the lower double-digit 
per-cent regime. 

Prestretching a material to the point where director reorientation sets in may allow to 
observe the electric-field-induced effect also in conventional non-swollen samples, as outlined 
theoretically [58]. Close to this critical point, the reorienting external stretching deformation 
on the one hand and the imprinted memory that anchors the director on the other hand 
approximately balance each other. Then a relatively weak additional external mechanical 
or electric held can result in a comparatively large response in the director orientation. 

Apart from that, other routes of regulation or actuation were demonstrated experimen¬ 
tally, for example photo- [59] 60] or temperature-induced [6T1 [64] deformations. Possible 
applications as soft actuators or even “artificial muscles” were discussed [65d69] . Large 
length changes of several 100 % were observed for loaded and unloaded samples under re- 
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Figure 6: Coupling between director reorientations and elastic deformations in nematic liquid crystal elas¬ 
tomers. The nematic director is indicated by the double-headed arrow, whereas the rectangular box il¬ 
lustrates the deformations of the elastomer. An initial ground state is depicted in the center. When an 
external stretch (marked by the vertical broad arrows on the left) is applied in the direction perpendicular 
to the initial director orientation, the coupling to the director can induce its reorientation into the stretch¬ 
ing direction. For practical purposes, the inverse effect would be beneficial, i.e. mechanical deformations 
induced by applying an electric field that reorients the director, as illustrated on the right. Unfortunately, 
the necessary field amplitudes turn out to be too high for practical purposes, except for materials swollen 
with conventional small-molecule liquid crystals. 


versible heating and cooling through the isotropic-nematic transition 53 [70]. For many 
applications, however, control parameters different from temperature are more desirable. 

On the theoretical side, we mention two principal routes of approach that reproduce 
the measured nonlinear stress-strain behavior as well as the detected deformations due to 
director reorientations. The first one is based on the concept of rubber elasticity [2]. It 
assumes Gaussian statistics for the step lengths along the polymer strands cm- As a starting 
point, it further assumes that the main role of the liquid crystalline building blocks is to 
introduce an anisotropy direction into these Gaussian statistics. Thus, instead of Eq. ([7]), 
the end-to-end vectors of the polymer strands between network points of the polymer mesh 
are assumed to follow the distorted distribution 

^(R) oc exp |-^R • (I-r a nn) ■ r| , (9) 

where n again indicates the nematic director, r a is an anisotropy parameter, I indicates 
the unity matrix, and L is an effective length connected to the actual contour lengths of 
the polymer strands. Open issues concerning this picture are how valid the assumption 
of the distorted Gaussian statistics generally is |[3j], and particularly how the orientational 
memory influences the statistics. A different microscopic approach was proposed in the form 
of randomly crosslinked anisotropically interacting rigid dimers [72], which reproduces the 
expressions for the theory of rubber elasticity on a more coarse-grained level. 
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Figure 7: Illustration of the concept of relative rotations in liquid crystalline elastomers. In the ground 
state (center) the director takes the orientation memorized by the architecture of the polymer network from 
the time of crosslinking. Relative rotations as a macroscopic variable describe how the director is rotated 
with respect to that memorized orientation. This can occur by a reorientation of the director relatively to 
the polymeric network (left), or a rotation of the network together with the imprinted memorized direction 
with respect to the actual director orientation (right). The variables of relative rotations test whether 
the results from both processes differ from each other and allow a formal inclusion in the corresponding 
energetic considerations. In all pictures, the strong-colored double-headed arrow indicates the actual director 
orientation, and the bright box symbolizes the polymeric network. On the left- and on the right-hand side, 
the light double-headed arrow represents the imprinted ground-state director orientation memorized by the 
network structure at the moment of crosslinking. 


The second route is macroscopic and based on symmetry arguments. It couples the 
continuum descriptions of liquid crystals and elastic materials ]73H53| . We noted above 
that a memory of the director orientation during synthesis is imprinted into the network 
structure of the polymer mesh. One central ingredient to many of the macroscopic studies 
is therefore the observation that rotations of the director out of its memorized equilibrium 
orientation cost energy. Such “relative rotations” between the director and the polymer 
network environment are included as additional variables (Z3E31IZBHZS]. First proposed for 
small deviations from the initial state, this concept was later generalized to the nonlinear 
regime mm- ft is schematically illustrated in Fig. [7| 

2.3. Ferrogels and magnetic elastomers 

Another example of externally tunable polymeric soft matter systems is given by mag¬ 
netic gels and elastomers [54]. Again, these materials combine the features of two different 
classes of complex fluids: those of ferro- and magnetorheological fluids [851 - 198] with those of 
crosslinked polymers [3]. 

The systems can be classified as composite magnetic hybrid materials [99H103] . They 
consist of superparamagnetic or ferromagnetic particles that are embedded in a crosslinked 
polymer matrix. Depending on the degree of swelling of such a magnetic elastomer with a 
solvent it may rather be called a magnetic gel [99, 104 , 105] . In the dry state, the materials 
can become hard and glass-like [1 99] , 
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Different sizes of the embedded magnetic particles, typically in the range of nano- to 
micrometers [84], can lead to qualitatively different material behavior. If the particles are 
smaller than ten to fifteen nanometers, the direction of their magnetic moment is not fixed 
with respect to the particle axes. Instead, thermal fluctuations can lead to a reorientation 
of the magnetic moment [106 ]. which is called Neel mechanism. For larger particles, reorien¬ 
tations of the magnetic moments mainly occur via the Brownian or Debye mechanism |M|, 
i.e. rigid rotations together with the whole particle. Still larger particles, from around a 
hundred nanometers, can feature multidomains of internal magnetization [1081 - 1112] . 

In the absence of an embedding polymer matrix, suspensions of such magnetic particles 
in a carrier liquid were stabilized in the form of ferro- and magnetorheological fluids [85H98] . 
Such fluids show properties of high practical relevance [ 113!. 1 114 ]. In particular, their dynamic 
flow behavior determined by their macroscopic viscosity can be tuned reversibly through an 
external magnetic field. This feature, commonly referred to as “magnetoviscous effect”, was 
observed experimentally and explained theoretically [871 SEE 11161 - 1123] . It is attributed on 
the one hand to a hindrance of particle rotations that would reorient the magnetization axes 
of the particles away from the aligning external magnetic field direction. On the other hand, 
it results from the formation of micro-aggregates in the presence of external magnetic fields. 
Typically, the second contribution, when present, exceeds the first one. The important point 
is that the adjustment of the viscosity can be achieved in a non-invasive way reversibly from 
outside. 

We turn back to the presence of an embedding polymer matrix in magnetic gels and 
elastomers. In analogy to the viscous behavior of magnetic fluids, for magnetic gels and 
elastomers the resistance to elastic deformations, quantified by the elastic moduli m, 
was experimentally probed. It turned out that the elastic moduli can be tuned reversibly 
from outside by external magnetic fields [84] 1125H13D] . This is one of the most outstanding 
properties of these materials. From an application point of view, it is for example interesting 
for the construction of novel externally tunable damping devices H3U or vibration absorbers 

m■ 

The analysis of theoretical minimum models demonstrated that the spatial distribution 
of the magnetic particles plays a central role and qualitatively determines the nature of this 
effect [13211134] . It was found that the tensile elastic modulus decreases or increases with 
increasing amplitude of an external magnetic field applied along the tension axis, depending 
on the particle arrangement. For instance, a regular simple-cubic or bcc lattice arrangement 
of the particle positions implied a decreasing tensile elastic modulus [ 1321 133]. In contrast 
to that, the modulus increased for an fee structure urn 

To a certain degree, the particle distribution can in fact be influenced during the synthesis 
of the materials. When a strong external magnetic field is applied before and during the 
crosslinking procedure, the magnetic particles were observed to form chain-like aggregates 
[841 113511139] . After the crosslinking has been completed, the field can be switched off with 
the anisotropic particle distribution persisting. The anisotropy has been locked into the 
materials. 

It has been demonstrated in experiments that this imprinted anisotropy can support 
the tunability of the elastic modulus by an external magnetic field pH]. In that case, the 
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largest effect of tunability was found when all three directions were aligned: the imprinted 
anisotropy direction, the additional external magnetic field to tune the mechanical prop¬ 
erties, as well as the direction of mechanical deformation to measure the elastic modulus 
[23]. This trend could be reproduced in a numerical study that assumed affine deformations 
of the material | 1140L 1141] . Although the degree of anisotropy in this case was significantly 
lower than that of chain-like aggregates, the effect was clearly observed. 

We remark that the assumption of affine deformations, which maps the macroscopic 
length changes linearly to all distances within the system, constitutes a major simplification. 
Yet it is often necessary to allow theoretical progress at all. Nevertheless, the more irregular 
the distributions of the magnetic particles within the samples become, the more this can 
lead to increasingly erroneous and even qualitatively incorrect results. Care has to be taken 
when this approximation is applied to realistic systems as has recently been pointed out on 
the basis of a simple dipole-spring model |142j . 

Allowing for non-affine deformations, an explicit numerical investigation on the impact 
of chain-like particle arrangements was performed using finite element simulations [ 143] . It 
showed that deviations from a straight chain architecture in the form of a zig-zag structure 
can qualitatively influence the behavior under external magnetic fields: the zig-zag angle 
determines whether the compressive elastic modulus decreases or increases with increasing 
amplitude of the external magnetic field. Even non-monotonic behavior is possible. 

The behavior within one linear chain was studied in a minimum dipole-spring model of 
hard spheres connected by elastic springs [144]. When the magnetic moment of the spheres is 
increased by an external magnetic held such that the hard spheres attract each other, they 
will come closer until they finally touch. Neglecting thermal fluctuations of the spheres, 
this process has the signature of a phase transition. It is of second order when the initial 
separation is small. If the initial separation is large, it is of first order, with the spheres 
clashing together in the final part of approach. A critical point separates both scenarios. 
Since the compressive elastic modulus along the chain diverges when the spheres touch each 
other, the process was called “hardening transition”. At the critical point, just before the 
spheres touch each other, the compressive elastic modulus along the chain vanishes. Also 
the behavior of free-standing hlamental chains was studied numerically [145 1 . [446] . In this 
case, various different states like compact, helicoidal, partially collapsed, simply closed, and 
extended open configurations were detected. 

A second fascinating property of magnetic elastomers and gels arises from their reversible 
deformations that can be induced by external magnetic fields. The spatial arrangement of 
the magnetic particles can qualitatively impact whether elongation or contraction occurs 
along the field direction, as pointed out by theory and simulations [ 147, 148] . An example 
is illustrated in Fig. [8] Apart from the magnetostrictive behavior [ 14711149TI153] . i.e. defor¬ 
mations induced by homogeneous external magnetic fields, particularly the “actuation” by 
inhomogeneous external magnetic fields [ 154 . 155] was investigated in experiments and by 
modeling. To minimize the free energy, the materials are drawn into regions of higher mag¬ 
netic field. This can lead to elastic deformations. After switching off the external field, the 
elastic energy is released and the probe can reversibly switch back to its initial state. It was 
demonstrated for a real sample that frequencies around 40 Hz can easily be followed [84 ]. 
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Figure 8: In magnetic gels and elastomers the spatial distribution of the magnetic particles can qualitatively 
influence the macroscopic behavior. Here, this is illustrated for the magnetostrictive behavior of a two- 
dimensional sheet in a homogeneous external magnetic field (see Ref. [147] V Left-hand side: for a random 
“gas-like” distribution of magnetic particles switching on the external magnetic field was observed to result 
in a contraction along the field direction [ I'47j . Right-hand side: in contrast to that, for randomly oriented 
doublets of magnetic particles, the same set-up was found to lead to an elongation along the field direction 

m- 


Consequently these materials again are ideal candidates for the construction of soft actuators 
[Ml 1156] or “artificial muscles” These applications involve dynamic properties, which 

so far have been addressed in only a few studies on the theoretical side [1581 - 1160] , requiring 
further investigation in the future. Apart from that, the use of magnetic gels as sensors 
indicating magnetic fields and held gradients by deformation was experimentally outlined 
[1610163] . Other studies discuss their application to control by external magnetic fields the 
amount and rate of drug release [ 1631. 163 ]. A further subject is the combat against cancer 
cells. For this purpose, magnetic particles are embedded into the corresponding cell tissue. 
A quickly alternating external magnetic held can be applied that continuously remagnetizes 
the particles. Through hysteretic losses during the remagnetization cycles, heat is generated. 
Such hyperthermic treatment [ 1121 1165H167] can destroy the cancer cells. 

For all these purposes, the degree of the so-called magneto-mechanical coupling should 
be maximized. A recent attempt in this direction was to not only embed the magnetic 
particles in the polymer matrix, but to directly crosslink them into the polymer network 
[104:. 1105 1 1168] . see Fig. [9] During sample preparation, this was achieved via surface- 
functionalization of the magnetic particles by attaching reactive groups on their surfaces. 
As a consequence, the magnetic particles themselves can serve as crosslinkers. Rotating the 
particles, the polymer chains attached to their surfaces get “rolled up” around the particles 
as investigated by a simulation study [169] . This “rolling-up” leads to restoring torques. An 
orientational memory of the particle orientations arises, implying an energetic penalty when 
the particles rotate relatively to their environment [144] , It can have a significant impact 
on the appearance of the system. For example, depending on the orientational memory, 
qualitatively different energetic ground states were obtained for a linear chain of ferromag¬ 
netic particles in a dipole-spring model [ 144] : ferromagnetic, with the magnetic moments 
aligned along the chain axis; antiferromagnetic-like, with neighboring dipoles rotated by 
180 degrees around the chain axis; or a spiral-like magnetization, with neighboring dipolar 
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Figure 9: Schematic illustration of two qualitatively different kinds of ferrogels implying different degrees of 
magneto-mechanical coupling, (a) In the classical materials the magnetic particles are simply embedded in 
the surrounding mesh of polymer chains that form the polymer matrix. For spherical particles, there may 
not be any permanent restoring torque induced by particle rotations, (b) In a new class of ferrogels, the 
polymer chains are covalently bound to the surface of the magnetic particles llOC 1105] . Thus the magnetic 
particles serve as crosslinkers. Rotations of the spherical particles inevitably induce deformations of the 
surrounding polymer mesh, leading to permanent restoring torques. In both cases, the density of polymer 
chains is significantly higher than indicated in the simplified schematics. 


moments rotated by an intermediate angle around the chain axis. Even stronger magneto- 
mechanical coupling might be achieved through the use of elongated instead of spherical 
magnetic particles [170] . 

2-4- Links between liquid crystalline and magnetic elastomers 

At the end of this section, we include a brief qualitative comparison between the two 
classes of materials considered above, i.e. liquid crystalline and magnetic elastomers and 
gels. We only concentrate on a few selected points that may be important when a mutual 
support of the two fields is considered. With the macroscopic theories for the two classes 
of material being partially mappable onto each other, also their observable behavior in 
analogous experiments may be so. 

Both materials combine the elastic behavior of polymeric substances with the anisotropy 
arising from an orientational ordering of their constituents. In the case of liquid crystalline 
elastomers, this anisotropy results from the average orientation of the liquid crystalline 
molecules along the director. For magnetic elastomers and gels, it is identified with the 
orientation of the magnetic moments in the equilibrium configuration. 

Depending on the procedure of synthesis, a more or less pronounced memory of the 
initial orientations at the crosslinking time is locked into the samples. Then relative rotations 
between this locked-in memory direction and the actual orientation cost energy. We outlined 
above the use of such relative rotations as a variable in the macroscopic characterization 
of liquid crystalline elastomers [SHI EH, ZU EX HZ], see Fig. [Tj and in a mesoscopic model 
of magnetic gels and elastomers [ 144] , Furthermore, relative rotations have already been 
considered in a macroscopic characterization of anisotropic magnetic gels [159] . A next step 
will be to connect the mesoscopic and the macroscopic descriptions of magnetic gels by an 
appropriate coarse-graining procedure. Similarities between the two classes of materials, i.e. 
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between liquid crystalline elastomers and magnetic gels, were pointed out on the macroscopic 
scale |171j . 

Generally, it is an interesting question, how the memory actually gets locked into the 
architecture on the molecular level. In other words, what structural changes occur for a 
polymer network in comparison to a conventional material without memory. The differences 
between the two classes of materials are of course large, already in view of the corresponding 
length scales. On the one hand, the liquid crystalline molecules are relatively small and can 
even be part of the polymer backbones in main-chain liquid crystalline elastomers. Therefore 
one would expect that an orientational memory should be reflected on the molecular level 
by the structure of the crosslinked polymer network. On the other hand, magnetic elas¬ 
tomers are real composite materials of mesoscopic colloidal magnetic particles embedded 
in a polymer matrix. Here, the spatial arrangement of the mesoscopic magnetic particles 
within the elastic matrix can be the source of the orientational memory. So the mechanism 
of memorizing may be quite different in the two cases. 

Nevertheless, based on the similarities in the macroscopic theories, it would be worthwhile 
to test the materials in those experiments in which the respective other material class features 
its outstanding properties. The experience acquired for one of the two material classes may 
help to guide the way also for the other material class. On the one hand, this means to 
clarify the question, whether an orientational coupling between strain deformations and the 
anisotropy direction also exists within anisotropic magnetic gels. In a next step, this would 
lead to the question, whether this property is likewise connected to a marked nonlinear 
stress-strain behavior as observed for liquid crystalline elastomers. On the other hand, in 
analogy to the elastic moduli that are reversibly tunable by external magnetic fields in the 
case of magnetic gels, a similar effect may be observable for liquid crystalline elastomers 
and gels. Here, rather an externally applied electric field would serve to tune the elastic 
properties. It would be satisfying to observe a mutual benefit between the two fields of 
liquid crystalline elastomers and magnetic gels. 

Moreover, instead of only tuning their properties, the materials could also be utilized 
as model systems for driven and active soft matter in genuinely non-equilibrium situations. 
We mentioned that both liquid crystalline elastomers and magnetic gels can be used as 
soft actuators by switching an appropriate external held [35], 09, 50] ES I5TH5711591 [601 62! 
EH] HU 11471114911157] . In a next step, these materials could be steadily driven by external 
fields so that they cannot reach an equilibrium state any more. Examples of driven soft 
matter will be addressed in the subsequent part of this review. In another step, they can be 
continuously “activated” and then be categorized as active soft materials. For instance, the 
potential of liquid crystalline elastomers for the construction of light-activated swimmers 
was experimentally outlined [ 172] . Similarly, an efficient production of light-activatable mi¬ 
croscopic artificial cilia was demonstrated via inkjet printing imii induced deformations 
of these slender artificial fibers [173] mimic the deformation cycle of natural cilia that are 
found on bacterial surfaces and serve for self-propulsion of the microorganisms [l74j. Like¬ 
wise, magnetically activated elastic filaments were used to construct propelling artificial 
microswimmers im We will come back to the migration of self-propelled particles and 
microswimmers in the later part of this review. Apart from that, on the theoretical side, 
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the macroscopic characterizations of gels featuring orientational degrees of freedom can be 
extended to include active stresses that result for example from chemical reactions. Such ap¬ 
proaches were used to identify and describe possible mechanisms for the motility of crawling 
biological cells |176H181] . 
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Figure 10: The flow profile corresponding to a planar linear shear flow is depicted on the left. It reveals 
the three different directions indicated on the right: the flow direction oriented along the flow velocity v 
(horizontal); the direction of the shear gradient V||v|| (vertical); and the vorticity direction along V x v, 
here oriented perpendicular to the shear plane. 


3. Driven soft matter 

Naturally, in each system thermal fluctuations lead to a sort of “internal drive”. For 
instance, colloidal particles are subject to Brownian motion, which adds to prevent them 
from sedimentation in a gravitational field and keeps them suspended [6]. Another example 
are experiments on colloidal particles thermally driven along the top of aligned membrane 
tubes that adhere to a substrate [182 1 1183] . An interesting crossover between exponential and 
Gaussian forms of the displacement statistics was observed in the latter context. Generally, 
when asymmetric ratchet-like conditions direct the internal thermal drive, net motion can 

arise mm- 

In contrast to that, we here focus on soft matter systems driven from outside. In our 
case this is achieved by externally imposed shear flows. The simplest situation is the one of 
an imposed planar linear shear profile as the one schematically indicated in Fig. 10 This 


geometry is characterized by three different directions: the flow direction, here horizontally 
oriented; the direction of the shear gradient, here vertically oriented; and the vorticity 
direction perpendicular to the shear plane, i.e. here pointing into the plane of the figure. 
In practice, the strategy to approximately realize this flow field in simple situations is to 
horizontally drive the top and/or bottom plate(s) of the confining cell with constant speed. 

We start by repeating simple effects that emerge for bulk-filling low-molecular-weight ne¬ 
matic liquid crystals. Increasing the complexity, we then briefly review aspects occurring for 
bulk-filling periodically modulated phases such as they appear in block copolymer systems. 
Finally we address localized structures in the form of vesicles. 

3.1. Flow alignment and director tumbling in nematic liquid crystals 

Similar to the previous chapter, we first turn to conventional low-molecular-weight liq¬ 
uid crystals deep in the nematic phase. Locally, the liquid crystalline molecules align on 
average along an ordering direction, the so-called director n [5], see also Sec. 2.1 In the 
hydrodynamic equations describing the macroscopic dynamic behavior of a nematic liquid 
crystal, the director orientation n couples to the hydrodynamic flow field v [5,186]. For an 


imposed flow field as the one depicted in Fig. 10 we now consider a bulk state of vanishing 
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Figure 11: Illustration of three dynamic modes of a nematic liquid crystal under an imposed planar shear. 
The nematic director is indicated by the double-headed arrow, whereas the orientation of the shear is given 
by the two arrows in the top right and bottom left in each panel, (a) In the case of “flow alignment” or “shear 
alignment” the director takes a steady orientation within the shear plane and with a finite tilt angle with 
respect to the flow direction, (b) When “wagging”, the director oscillates within the shear plane around an 
orientation of usually non-vanishing tilt with respect to the flow direction, (c) During “tumbling” motions 
the director makes full rotations in the shear plane. 


director distortions Vn = 0. Then the hydrodynamic equation for the director orientation 
becomes [186, IT87 j 

r)r\ 

n + A [A ■ n — (n • A • n) n]. (10) 


— = n 

dt 


In general, the two gradient tensors f2 = [(Vv) r — (Vv)]/2 and A = [(Vv) T + (Vv)]/2 
contain the rotational and elongational effect of the flow field v, respectively. Here Vv again 
denotes a dyadic product and the superscript T its transpose. A convective contribution 
v- Vn vanishes in our case due to the assumed spatially homogeneous director orientation 

Vn = 0. 

We can see from Eq. (10) that, besides the shear rate set by the imposed flow field, 
only one parameter A determines the dynamic behavior. It can be written as the ratio 
between two viscosities and thus represents material properties [5] H"87| . Depending on the 
magnitude of this parameter, the director can show qualitatively different types of behavior 
in the shear flow. On the one hand, for | A| > 1, the director reaches a steady-state orientation 
within the shear plane PUHSM. This situation is typically called “flow alignment” or 
“shear alignment 


a). The orientation angle with respect to the 
flow alignment angle”, remains constant and is determined by the 


and is depicted in Fig. 11 
flow direction, i.e. the 

parameter A. On the other hand, the director does not find a steady-state orientation, but 
continuously rotates in the shear plane, for |A| < 1. This dynamic state is usually referred 
to as “tumbling” HBsmsa, see Fig. [Tl] (c), and consequently A is often referred to as the 
tumbling parameter. Using the Fokker-Planck approach, these two dynamic regimes were 
connected to properties on the molecular level |191H194] . The above considerations imply 
that the shear rate of the imposed flow field is low enough so that the degree of nematic 
ordering, see Eq. ([Tj), is not significantly altered from its equilibrium value. 

Director tumbling was also observed and discussed for mixtures of liquid crystals ra 
and for liquid crystalline polymers |1871 flUS] . For certain parameter values, a transition 
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Figure 12: Schematic illustration of the two-dimensional patterns shown by micro-phase separated linear 
diblock copolymers. The temperature must be low enough so that micro-phase separation can take place. 
In the series from left to right, the fraction of one kind of polymer on the diblock copolymer chain with 
respect to the other one increases. This is depicted by a typical chain in each panel. If one species strongly 
dominates by its volume fraction, the state is disordered, i.e. not micro-phase separated, as in the most 
left and right examples. When the volume fractions become more balanced but still one species dominates, 
hexagonal textures can be observed as in the second examples from left and right. In three dimensions these 
would correspond to hexagonally arranged cylinders elongated perpendicular to the plane of the picture. 
Lamellar structures (center) appear for approximately equal volume fractions of the two blocks and are 
likewise obtained in three dimensions. 


was found with increasing shear rate from tumbling to flow alignment via an intermediate 
third dynamic state called “wagging” [1871 1197] , In this third state, which is included 
in Fig. 11 (b), the director oscillates back and forth instead of performing full tumbling 
rotations. Still further dynamic modes were identified later, including chaotic types of 
motion [198] . 


3.2. Alignment and reorientations in periodically modulated phases 

At our next level of complexity, we concentrate on materials that in their ground states 
cannot be considered as spatially homogeneous any more. Instead, they feature regular 
periodic modulations in their density or in the concentrations of their constituents. An 
illustrative example are block copolymer melts or solutions. 

To understand the situation in block copolymer melts, it is instructive to recall the 
problem of mixing of two polymers of a different kind. This scenario is addressed by the 
Flory-Huggins theory [3]. For low-molecular-weight liquids, a spontaneous mixing is usually 
driven by a gain in translational entropy: each molecule is now provided with the extended 
volume of the whole mixture, not only with the one of its initial single component liquid. In 
polymers, however, this effect is strongly reduced. The monomers are chemically trapped in 
the polymer chains. Only each polymer chain as a whole can increase its translational en¬ 
tropy. Therefore the gain in translational entropy under mixing is significantly lower than for 
low-molecular-weight liquids. Instead, the inter-species interactions between the monomers 
of different kinds in comparison to the intra-species interactions between monomers of the 
same kind determine whether the polymers mix or form two separated phases. In most 
cases, this competition favors the separated state. Most pairs of polymers do not form 
homogeneous mixtures at ambient temperatures [3j. 

The situation becomes interesting in block copolymers. Here, the different polymer 
chains that tend to demix under the previous conditions are combined into one molecule. 
For example, a diblock copolymer chain is composed of two chains of chemically different 
polymers that are covalently bound to each other at one of their ends. Due to the covalent 
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Figure 13: Illustration of the three principal orientations of a lamellar sheet-like structure under an imposed 
planar shear: (a) parallel, (b) transverse, and (c) perpendicular. 


bonds between the blocks, a macroscopic demixing is not possible for block copolymers. 
Consequently, demixing can only occur on the length scale of a single block in the form of 
a micro-phase separation. These demixed blocks can then arrange into spatially periodic 
patterns. Depending mostly on temperature and the volume fractions of the two different 
blocks, lamellar phases, hexagonally arranged cylinders, bcc or fee (in solutions) structures, 
and more complicated textures are obtained 0 US H99H208]. For approximately equal 
volume fractions of the two blocks, a lamellar micro-phase separated state emerges. It 
features a layered arrangement of the different blocks and was depicted in Fig. [l] Patterns 
expected for two-dimensional arrangements of linear diblock copolymers are depicted in 


Fig. 12 


When such materials are exposed to shear, an orientation of the periodic structures with 
respect to the shear velocity and its gradient has been observed in numerous experimental 
studies (204 , 20511209H228] : see Fig. 10 for the shear geometry. The most illustrative example 
is the one of lamellar layered textures [2U41 120511209H222L1224112251122711228] . For a simple 
planar linear shear flow, three principal orientations of the lamellae are possible as illustrated 
in Fig. 13: a “parallel” orientation with the lamellae normals parallel to the direction of the 
shear gradient; a “transverse” orientation with the lamellae normals parallel to the direction 
of the shear flow; and a “perpendicular” orientation with the lamellae normals perpendicular 
to the two previous cases, i.e. parallel to the vorticity direction. 

Interestingly, when a steady linear shear flow was applied to a lamellar block copolymer 
sample, the signature of the transverse orientation was observed at low shear rates |222j . 
This is surprising because it leads to a significant distortion of the lamellar layers. Switching 
to higher shear rates, a reorientation process occurred with the layer normals now lying in 
the plane spanned by the directions of the shear gradient and the vorticity |222j. This 


implies a combination between the parallel and the perpendicular alignment, which was also 
observed in other studies [2111 212j. However, the steady shear could not induce a perfect 
reorientation. Even more, the steady shear led to defects in previously well-aligned lamellar 
samples for example in the form of focal conical textures [210] , It turned out that large- 
amplitude oscillatory shear was much more effective to align these samples than steady shear 
flow. Thus the orientational effects of oscillatory shear were studied extensively. 
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plate-plate cone-plate 

Figure 14: Schematic illustration of two principal types of shear rheometer, (a) In the plate-plate geometry, 
the sample is sandwiched between two parallel plates, one of them rotating relatively with respect to the 
other around their common axis. The displacements on the plate surface and thus the shear amplitude 
increase linearly from the rotation axis in the center to the outside, (b) A cone-plate geometry counteracts 
this effect because the sample thickness likewise increases linearly from the central rotation axis to the 
outside. In both cases, steady shear (darker single-headed arrow) or oscillatory shear (brighter double¬ 
headed arrow) can be applied. The arrows also coincide with the local direction of the shear velocity, see 
Fig. [10} 


For large-amplitude oscillatory shear, reorientation transitions in lamellar phases were 
observed by varying the shear frequency or the shear amplitude [2051 2091 212112211 12241 
12251 1228 ]. Even a “double-flip” reorientation process can be found for increasing shear 
frequency. At low frequencies, a parallel orientation is often obtained [211112T21 2T4 . 2161 - 
12191 22 1 . 225 J. An increase in the shear frequency can lead to a flip to the perpendicular 
orientation [214, 216-219, 224. 225J. Very high frequencies, however, can induce a second flip 
back into the parallel geometry [2T3112141 218 . 21911224] . Despite the significant distortion 
of the lamellae, the transverse orientation was also observed in the case of oscillatory shear 

[mi Ena E2i]. 


Experimentally, the investigations can be performed in a Couette cell or using rheometric 

In the plate-plate geometry, the sample is sandwiched 


devices as displayed in Fig. 14 


between two parallel plates that rotate relatively with respect to each other around their 
common axis. At the surface of the rotating plate, the displacement increases linearly 
with increasing distance from the rotation axis in the center. Thus effects of different 
shear amplitude, defined here as displacement versus sample thickness, can be investigated 
simultaneously. A cone-plate geometry counteracts this effect. Here, one plate is replaced by 
a cone, leading to a linearly increasing sample thickness from the rotation axis in the center to 
the outside. Generally, the principal lamellar orientation in the sample cell can for instance 
be determined by small-angle x-ray or neutron scattering [2041 205112091 - 122511227 , 228]. From 
the results of different beam directions, conclusions can be drawn about the orientational 
state. 

Understanding the underlying reasons for the orientational behavior under shear is a 
complicated task and also material dependent. Several different effects connected to various 
length scales play together and determine the overall behavior. For example, the boundary 
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of the sample cavity can have an aligning effect on the structure [229] ; the micro-phase 
separated textures are distorted by shear deformations, leading to elastic strains 
structural fluctuations can couple to shear deformations [ 233] ; viscosity contrasts are present 
in each sample due to the stacking of the molecules and due to the material properties of the 
different blocks [2131 [214] ; the molecules are elongated along the layer normal, which sets a 
preferred path of molecular motion in particular for entangled polymers at higher molecular 
weights 123 ; shear can induce a tilting of the elongated chains away from the layer normal, 
leading to a restoring torque [2201.123411236] : and as we have seen in the previous section on 
low-molecular-weight liquid crystals, elongated molecules by themselves can already show 
various types of dynamic alignment behavior. Apparently, different frequencies and ampli¬ 
tudes of the applied shear can be used to preferably address some of these ingredients. In 
this way, their relative impact on the collective dynamic response can be modified and the 
overall appearance of the sample can be tuned. 

Other materials like common smectic (layered) low-molecular-weight liquid crystals [ 2331 
I237H239] . see Fig.g(c), or lamellar phases of surfactant solutions [233112401 - 1244] can similarly 
feature orientational effects due to their “sheet-like” structure. Reorientations under shear 
can also be observed in structural phases different from the lamellar one [2M 12U51 m 
12271 1245 . 246] . Besides, for hexagonally arranged cylinders as well as for spheres ordered 
in bcc- and fcc-textures, a sliding of whole structural sheets over each other was suggested 
[2(751 1245112TW252] . 

Applying shear to block copolymer systems is an issue of high practical relevance. Due to 
the reorientation effects, it can serve to heal multi-domain textures and obtain a monodomain 
structure [2T011221. 223 . 233 .. 125TI1253] . Apart from shear, alignments can also be induced 
by other routes, for example by applying external electrical fields to materials that show a 
dielectric contrast between their blocks [22911254H267] . From a technological point of view, 
block copolymers possess a high potential for applications in the nano-sciences [26811270] : 
they provide regular periodic arrays of tunable material properties on the molecular length 
scales of the different blocks. This can be exploited for nano-lithographic processes, to 
produce regular arrays of nano-dots or nano-wires, or to synthesize nanoporous materials. 
Monodomain structures are often beneficial for this purpose. Therefore understanding and 
clarifying the possible routes to generate them is a task of central importance. 


3.3. Vesicles in shear flow 

In the above cases, we considered the properties of bulk-filling phases under shear flow. 
This is different from the following examples, where we focus our attention on localized 
finitely-sized objects such as closed membranes and vesicles. In contrast to a bulk-filling 
periodic phase, complete rotations of a limited object in a liquid environment is easily pos¬ 
sible. Together with the deformability of the considered entities, this leads to new dynamic 
states. 

The systems that we focus our attention on are based on amphiphilic molecules in an 
aqueous environment [ 271 . [272], Such molecules typically feature a hydrophilic (“water- 
loving” ) head group and one or more hydrophobic (“water-fearing”) tails, usually based on 
hydrocarbon chains; see Figs, g and 15 (a). 
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Figure 15: Schematic examples of amphiphilic (surfactant) molecules and illustration of some of the ag¬ 
gregates that they can form in an aqueous solution, (a) Two extreme examples of typical structures of 
amphiphilic molecules are displayed. Both have a hydrophilic (“water-loving”) head group, typically in 
the form of some polar or charged molecular groups. The tails are usually based on hydrocarbon chains 
and hydrophobic (“water-fearing”). In the first example, the head group is relatively bulky and there is 
only one shorter tail. In the second example, there are two tails of approximately the same diameter as 
the head group. The formation of aggregates in an aqueous environment in both cases is possible above a 
certain concentration. It leads to a screening of the hydrophobic tails from the aqueous environment by the 
hydrophilic head groups, (b) Due to their geometric properties, the molecules of the first kind favor the 
formation of micellar aggregates of high surface curvature, (c) The molecules of the second kind, due to 
their geometric shape, are suitable to generate flat bilayer membranes. Thermal fluctuations are of course 
present but not indicated here. To close the membranes at their outer ends, some arrangements of high 
surface curvature and thus higher energy are necessary, (d) For large membranes, this can be avoided by 
forming closed surfaces in the form of vesicles. The ratio of their radius versus the membrane thickness in 
reality can be much larger than depicted here, so that locally the surface curvature of the vesicle becomes 
negligible. See also the indicated sizes in Fig. [T] 


To reduce the contact with the surrounding water molecules, the hydrophobic chains 
above a critical concentration tend to form intermolecular aggregates. This is again a sort of 
micro-phase separation. A simple example is depicted in Fig. [l5](b) in the form of a micelle. 
In this aggregate, the hydrophobic chains can hide away from the water molecules and form 
the core of the micelle. They are screened from the water molecules by the outward-pointing 
hydrophilic head groups that search the contact to the aqueous environment. 

Both, the concentration and the geometry of the single amphiphilic molecules determine 
which sort of intermolecular aggregates form. If the head groups are very bulky and if 
there is only one short chain per molecule, the formation of spherical micelles as depicted 
in Fig. 15 (b) is supported. They have a high curvature of their surface. In contrast to 
that, if the hydrophilic head and the hydrophobic tail groups have approximately the same 
diameter, the formation of aggregates with low curvature of the surface is beneficial. This 
situation can for example be observed for bilayer membranes. In this case, the amphiphilic 
molecules are arranged in two adjacent layers stacked on each other as illustrated in Figs, [l] 
and 15 (c). Here, the hydrophobic chains point to the inside of the two stacked layers and 
expel the water molecules. The hydrophilic heads point to the outside and screen the inner 
region from the aqueous environment. 

At its lateral rims, the bilayer membrane would either have to expose some of the hy¬ 
drophobic chains to the aqueous environment or show an unnaturally high curvature of its 
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Figure 16: A schematic vesicle is half-cut to show the membrane location. Within the plane of the membrane, 
a shear stress is applied to the membrane as indicated by the two arrows. If the membrane does not sustain an 
in-plane shear stress but behaves like a two-dimensional liquid, the vesicle is called a “fluid vesicle”. Likewise, 
it is referred to as a “viscous vesicle” when the effect of the intra-membrane viscosity is emphasized. In 
contrast to that, if the membrane does sustain in-plane shear stresses, the vesicle is called an “elastic vesicle”. 


surface. To avoid such high-energy regions at the lateral rims, the bilayer can bend on length 
scales large compared to its thickness and form a closed surface. The resulting closed bilayer 
membrane is called a vesicle and also indicated in Figs. [I] and 15 (d). Typically, an aqueous 
solution is found on both sides of the bilayer membrane, i.e. outside the vesicle and also 
enclosed on the inside. 

Interesting dynamic states arise due to the typical constraints following from the archi¬ 
tecture of vesicles. On experimental time scales, their inner volume is generally conserved 
[ 27211273 ] , Most importantly, also the surface area of vesicles remains practically constant 
[ 27211273 ] . This is in contrast to simple liquid drops and qualitatively affects their shape 
and dynamics [ 274 . 275 ]. Apart from that, a possible viscosity contrast between the inner 
encapsulated fluid and the outer fluid impacts the dynamic behavior [ 274 , 127614302] . The 
same is true for a possible friction between the two layers of amphiphilic molecules that form 
the bilayer membrane and may slide over each other during certain dynamic modes [ 272 ] . 

On the one hand, the membrane may correspond to a two-dimensional liquid that does 
not sustain in-plane shear stresses, see Fig. 16 Then the vesicles are called “fluid vesicles” 


[271'. [272] . We concentrate on the dynamics of such objects in simple linear planar shear 
flow as introduced in Fig. [10] If there are no viscosity contrasts between the liquid on the 
inside and on the outside of the vesicle, a dynamic state called “tank-treading” is forecast 
by modeling and observed experimentally [27514278112811128311285H287112911129211297112981 
1303. 304j . It is depicted in Fig. 17 (a): due to the elongational part of the shear flow the 


vesicle becomes oriented within the shear plane, but with a constant inclination angle to 
the flow direction; the rotational part of the shear flow induces a rotational motion of the 
membrane around the vesicle interior. Increasing the viscosity contrast at not too high 
shear rates, a bifurcation to a different type of motion called “tumbling” is predicted from 
modeling and observed experimentally [275H279112821128311285H287L1291112921129711298j . It 
is illustrated in Fig. 17 (c). Here, the vesicle as a whole rotates in the shear flow, instead 
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tank-treading swinging tumbling 


Figure 17: Illustration of three dynamic modes of a vesicle under an imposed planar shear flow. The 
vesicle is indicated by the ellipsoid, (a) While “tank-treading”, the vesicle takes a steady finitely inclined 
orientation with respect to the flow direction within the shear plane. Only the vesicle membrane performs 
a continuous rotational motion along its contour, (b) During “swinging” motions (also “trembling” or 
“vacillating breathing”) the long axis of the vesicle oscillates within the shear plane around an orientation 
that is usually finitely inclined with respect to the flow direction, (c) When “tumbling” the vesicle performs 
full rotations in the shear plane. 


of only its membrane tank-treading around its interior. Thermal fluctuations can induce 
tumbling already for vanishing viscosity contrast at low shear rates [305j. At higher shear 
rates, “swinging” (also “trembling” or “vacillating breathing”) was observed instead of pure 
tumbling at least for intermediate viscosity contrasts |2751128211284H28711291112921129711298] . 
In this state, the long axis of the vesicle does not perform full rotations as in the tumbling 
mode, but only oscillates up and down as indicated in Fig. 17 (b). Summarizing, a tendency 
towards the right-hand side in the series of panels (a)-(c) in Fig. 17 is generally supported 
by an increasing viscosity contrast between the outside and the inside of the vesicle and a 
decreasing shear rate. Yet, the situation can be more involved in each individual case. More 
complex dynamic modes were identified recently [ 297 , 299j, and the qualitative influence 
that thermal fluctuations can have were pointed out PSHUSMDE]. 

Apart from that, the impact of an explicitly viscous membrane was investigated [285 . 
130911311] . These “viscous vesicles” generally show a reduced inclination angle during tank- 
treading. Increasing the membrane viscosity can further induce a transition from tank- 
treading to tumbling or swinging. 

On the other hand, the membrane may for example be polymerized so that it does sustain 
in-plane shear stresses pm 1272] . see Fig. [L6] Such objects are referred to as “clastic vesicles” 


[ 312] , The basic features of the dynamics of fluid vesicles in linear shear flow are recovered 
both from modeling and experiments: a tumbling motion at low shear rate and a tank- 
treading motion at high shear rate emerge [290112931, 305 , 131311315] . However, during the 
tank-treading motion, the orientation angle with respect to the flow direction oscillates due to 
the shape-memory [ 2901 293 . 3051131311315] . This is referred to as “swinging”. Furthermore, 
if the viscosity inside becomes large compared to the outside viscosity, an intermittent type 
of motion is observed as a combination between the swinging tank-treading and the tumbling 
dynamics [2901 12931 [3051 SEE 1313] . 

From a biological point of view, elastic vesicles form an important topic because red 
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Figure 18: Schematic illustration of different possible shapes adopted by vesicles in a Poiseuille flow. The 
snapshots have to be considered as two-dimensional cross-sectional views through the center of a tube-like 
channel, (a) Profile of the Poiseuille flow across the channel, (b) Elongation of a discoidal vesicle along the 
fluid flow, (c) Parachute-like shape and orientation, (d) Slipper-like configuration. Gray bars indicate the 
channel walls. 


blood cells fall into this category [31311316j . To understand the transport of red blood cells 
in blood vessels, the dynamics of elastic vesicles subjected to Poiseuille flows in cylindrical 
tubes, see Fig. 18 (a), is studied extensively. Whereas fluid vesicles of initially discoidal shape 
typically elongate in the Poiseuille flow, see Fig. 18 (b), or show axisymmetric “parachute” 
shapes in strong flow fields, see Fig. 18 (c), elastic vesicles usually adopt a parachute-like 
deformation or an asymmetric “slipper”-like shape as indicated in Fig. 18 (d) [3T2, 317] . 


Many vesicles together in a tube were found to adopt disordered distributions with discoidal 
shapes, ordered states of aligned parachutes, as well as slippers arranged in a zigzag-like 
fashion [31811320] depending on the flow velocity and packing fraction. Hydrodynamic lift 
forces keep the vesicles away from the walls [ 291 . 132114326] . In addition to that, further 
dynamic effects have been discovered, for example a wrinkling deformation of the elastic 
membrane in shear flows [M132711525] . 

Analytical approaches included the above-mentioned constraints following from the vesi¬ 
cle architecture and were able to analyze the resulting dynamic properties [27411275112791.12841 - 
129011293112931 [295112961 [3001 [3CT5113081132T1 [3281 [33T] . Due to the curved vesicle surfaces, 
elements from differential geometry are necessary to perform these calculations. Apart from 
the analytical approaches, two simulation routes were particularly successful to investigate 
the dynamic behavior of vesicles in a liquid environment. One of them is particle-based and 
referred to as multi-particle collision dynamics [289112911130211309H31211318113201133211333] . 
It includes the surrounding and encapsulated fluid particles, modeling the many-particle 
interactions in an effective way. Still, it correctly obeys the physical conservation laws. The 
second route is a phase held model [276H27811280112941133411335] . It is based on an aux¬ 
iliary continuous scalar order parameter held. Two different constant values of this scalar 
order parameter are used to distinguish between the inside and the outside of the vesicle, 
for example —1 and +1, respectively. At the location of the membrane, the order parameter 
held rapidly but continuously changes between these two values. In this way, the boundary 
of the vesicle can easily be identified. It can be effectively tracked by simply advecting the 
order parameter held with the hydrodynamic how of the system. 
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A third route to the problem was outlined in the form of a density-field approach [336]. 
For periodic bulk-filling textures, such a procedure had for example been introduced for 
block copolymer structures [201 , 1230H232112471124811337H341] or for crystalline materials via 
the phase field crystal approach |342l 1349] . These kinds of formalism are based on a single 
scalar real order parameter naturally given by the density field. They apply on diffusive time 
scales. A corresponding free energy density to obtain localized nonperiodic structures was 
introduced in a phenomenological way and the resulting density-field description was coupled 
to hydrodynamic flow fields [ 336] , In this way, in two spatial dimensions, characteristics of 
tank-treading were observed in linear shear flow, and folding into the parachute form as well 
as an elongation in a Poiseuille channel flow were obtained [336]. A challenging task for 
the future is to derive such density-field approaches starting from a more microscopic and 
particle-resolved basis. 


31 
















4. Active soft matter 


So far we have been concerned with “passive” systems that at most were driven from 
outside. We now turn to “active” or “self-driven” materials. It is not always straightforward 
to clearly distinguish between such “active” soft matter [35UU354] and “externally driven” 
soft matter as considered in the previous section. The term “active” generally implies that at 
least some components of the system feature an individual “internal drive”. Unfortunately, 
this definition requires a sort of restricted point of view. There is always a residual coupling 
of the internal driving mechanism to the outside world. We understand this statement by 
considering examples of active self-propelled particles. 

One case of particularly well studied artificial self-propelled microswimmers is given by 
dispersed colloidal Janus particles |355fl361j . In general, colloids are made of mesoscopic 
particles or droplets of sizes between 1 nm and 10 — 100 /ini, see also Fig. [lj dispersed in 
another substance [SHE]. Everyday examples are ink, paint, and milk. To stabilize colloidal 
suspensions, coagulation of the particles must be hindered. For this purpose, the refractive 
index of the colloidal particles and their environment should be matched to each other to 
reduce the van-der-Waals interaction. Furthermore, electric charges on the particles can lead 
to stabilization due to repulsion. Apart from that, also steric stabilization is possible, when 
the colloidal particles are covered by polymer brushes: close contact between the particles 
would restrict the possible configurations of the polymer chains, which induces an effective 
entropic short-range repulsion. Micron-sized colloidal particles offer the advantage that they 
can easily be tracked by optical methods. 

There are different ways to cover one side of the surface of a colloidal particle by a 
material of significantly different properties than the other side [T6HTK1 1362 ]. Such colloidal 
particles are then called Janus particles. This can be exploited, e.g., to selectively heat 
the stronger absorbing side by laser light irradiated from outside ES7J. As a result, a self- 
propelled motion can be induced via thermophoretic effects: the particle via its asymmetric 
light absorption generates a temperature gradient around itself; this leads to hydrodynamic 
stress gradients that can set the surrounding fluid into motion, which effectively leads to 
particle motion [363] . A further possibility is to enforce demixing of a binary fluid on the 
heated side [ 358, 359.1361] or to selectively catalyze chemical reactions only on the surface 
of one side of the particles [ 356, [360] . The arising concentration gradients can lead to 
self-propulsion via diffusiophoretic effects analogous to the case of thermophoresis [363]. 
However, the fuel to drive these chemical reactions needs to be provided from outside, as 
must be the irradiated laser light in the case of thermophoretic motion. 

Granular hoppers on a vibrating plate form another model system to study self-propelled 
motions [363H567] . Here, typically the motion in the plane parallel to the surface of the plate 
is observed. However, this motion is initiated and driven from outside by the externally 
tuned vibration. Finally, even such units as swimming bacteria [36811371] . crawling cells and 
amoebae [37211374] . or molecular motors that can set filaments into motion [37511376] are 
not completely independent of their surroundings. At some point they need to take up the 
food or fuel provided by their environment. 

We thus further narrow our present definition of active motion. The internal drive of the 
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Figure 19: Illustration of the Vicsek model 1377< in two spatial dimensions. All particles, marked by the 
dark bullets, self-propel with the same speed vq. This is indicated by the velocity vectors attached to the 
bullets that all have the same magnitude. At each discrete numerical time t , the velocity orientation dgf) 
of each particle i is updated: it is set to the average velocity orientation of all particles located within a 
distance tq from the ith particle; furthermore some orientational noise is added. During each time step At, 
the particles move along their velocity orientations by a discrete distance v^At, here shown for a unit time 
step At = 1. 


active units may be induced from outside. However, their motion must break a symmetry 
of the system. In other words, there must be a degree of freedom in the direction of induced 
motion that is not prescribed by the external driving. For example, a vertical vibration of 
the substrate can induce the motion of granular hoppers on a horizontal plate. However, 
the direction of motion within the plane of the plate surface is not specified. In that sense, 
the active particles are free to choose their migration direction. 

In the following, we are mainly interested in collective effects that arise when many active 
and self-propelled particles act together. Several aspects of this complex behavior seem to 
have general character and can already be understood from minimal model systems. There¬ 
fore, we increase the level of complexity of the systems under investigation. We start from 
the simplest case of point-like self-propelled particles featuring a constant self-propulsion 
velocity. Although this is an idealized situation, basic principles of the collective behavior 
are revealed. In a next step, we include systems where steric interactions occur, e.g. dry 
granular systems. After that, we outline the consequences of replacing the constant self¬ 
propulsion velocity by an active driving force. Artificial microswimmers like self-propelling 
colloidal Janus particles and biological swimmers in the form of bacteria are addressed. The 
influence of hydrodynamic interactions on their behavior is discussed. We then outline the 
role of particle defornrability. Finally, we give examples of studies dedicated to the collective 
behavior of animals. 


4-1. Point-like self-propelled particles 

The most prominent and fundamental model introduced in the context of idealized point- 
like particles is the one by Vicsek et al. [ 377 ]. which is illustrated in Fig. 19 and explained in 
the following. In this model, N particles self-propel in a two-dimensional plane of periodic 
boundary conditions. As a key ingredient and major simplification, the magnitude of the 
individual velocities of all particles, i.e. their speed, is assumed to be equal to a constant 
v 0 . It remains identical and unchanged for all times. The current orientation of the velocity 
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Figure 20: Schematic illustration of the order-disorder transition in the Vicsek model and its variants. The 
transition can be induced by increasing the characteristic noise amplitude or by decreasing the mean particle 
density. In the ordered state (bottom left inset) the particles on average migrate collectively into a common 
direction, whereas in the disordered state (bottom right inset) coherent particle motion does not occur. 
The order parameter P is calculated from the magnitude of the sample-averaged velocity orientations, see 
Eq. (131. Close to the transition, spatial inhomogeneities are usually observed in the form of density bands 
that travel perpendicularly to their elongation direction through a diluted disordered background (top right 
inset). 


vector of the ?’th particle can be parameterized by an angle i = 1,..., N. This orientational 
angle is adjusted at each time step to the mean of the velocity orientations of all particles 
j that are located within a spherical environment of radius ro around the zth particle: 


$i(t + At) = {fij{t)) r o + 


( 11 ) 


Here, (...) ro describes this average over all particles within the sphere of radius ro around 
the ith particle, t marks the time and At sets the discrete time step. Furthermore, the 
angular noise rji(t) is originally taken as a random number with uniform probability out of 
a centered interval at each time step. After each time step At, the particle positions r* are 
updated as 


r i{t + At) = r i(t) + Atv 0 


cos[#i(t)] \ 
sin [$i(t)] J ' 


( 12 ) 


As a result of the competition between the alignment with the local environment and 
the stochastic noise, a phase transition is observed. This is a transition between a phase of 
disordered motion of zero net particle flux on the one hand and a phase of ordered collective 
motion on the other hand. In the second case, all self-propelled particles migrate on average 
collectively into the same direction. The two phases are schematically indicated in the 


bottom insets of Fig. 20 


An order parameter to quantify this order-disorder transition is given by the magnitude 
of the sample-averaged velocity orientations, 


P(t) = 


]_ ( COS [&i(t)] \ 

N sin [i?j(t)] J 
1=1 
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where the velocity orientations are parameterized by the orientation angles dj (i — 1, N ) 
[see Eqs. and ( fl2| ) as well as Fig. [T9] . The magnitude of the order parameter is P = 1 
in a completely ordered state of all particles collectively migrating into the same direction, 
and P = 0 in the disordered state. As indicated in Fig. [20} the phase transition from the 
ordered to the disordered state of motion can be induced by decreasing the mean particle 
density or increasing the characteristic noise amplitude. 

This transition was studied in detail over the past few years [37811380] . In simulations of 
large systems, it was found that the transition is of first order [37811381] . The discontinuity 
in the transition is apparently related to spatial inhomogeneities that arise in the particle 
density around the transition point. Density bands emerge that tend to collectively migrate 
perpendicularly to their elongation [3791 380 . 13821 1383 j as indicated by the top inset in 
Fig. [20} Due to this directed collective migration, they can pick up further particles from 
the environment of disordered motion. Such a process further increases the density within 
the band. In this sense, a kind of self-supporting mechanism develops [384] . Also in real 
experiments, traveling density bands have been observed [375] . 

The nature of these traveling density bands was discussed in the framework of solitons 
[3791 [380] , Typically the density bands feature a sharp front and an extended tail. Recently, 
their behavior under head-on collisions is increasingly investigated for point- and non-point- 
likc particles [38511.387] . Penetration of colliding density bands and recovery after collision 
have been observed [385 . 1386 ]. It has been demonstrated that they are obtained as different 
propagating solutions of conventional continuum models for self-propelled particle crowds 
in the form of multiple parallel density bands, single solitary bands, and single 


active droplets 

Several variants of the Vicsek model were pointed out and analyzed [ 39lj . For example, 
the effects of metric-free alignment interactions [392] were discussed. In three spatial dimen¬ 
sions, the traveling density bands were recovered in the form of migrating density planes 
[379] , Other studies replaced the discrete nature of the dynamic equations Eqs. 0 and 


(12) by differential equations. In particular, the discrete averaging process (...) ro in Eq. 0 
was replaced by a more continuous functional form [3931 - 1395] : 


d'&ijt) 

dt 


dU dri 

■ w + nt), ^ = *.„ 


COS [$j(t)] 

sin [#;(£)] 


i = 1, ■, N. 


(14) 


For simplicity, the orientational noise Tj(f) is assumed to result from a Gaussian white 
process. The continuous function U is based on pairwise alignment interactions between the 
particles. For example, the functional form 


N 


U(Ti,...,Tjr,0i,.. .,&n) = - V Q t r 0 - ll r i - Tj ||) COs(l? 4 - $,) 


(15) 


i,j =1 
i<j 


again leads to pairwise velocity alignment for particles closer to each other than the distance 
r 0 . Here, 0 represents the Heaviside step function. The snapshot in the top inset of Fig. 20 
was obtained from a numerical calculation following a procedure along these lines. 
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Figure 21: Schematic illustration of polar alignment rules on the left and apolar alignment rules on the 
right. Self-propelled particles are indicated by rectangular boxes, while the arrows inside the boxes mark 
the direction of active drive. For polar alignment rules, the particles always tend to align their velocity 
orientations into the same direction. In the case of apolar alignment rules, which may be induced, for 
example, by steric interactions of elongated particles, an antiparallel orientation of the velocity vectors is 
equally preferred. 


Using such an approach, the situation of a binary mixture was considered for different 
rules of the inter-species velocity alignment [395j. Starting from the particle picture, con¬ 
tinuum equations for the one-particle probability densities were derived [395] within the 
Fokker-Planck framework [ 191 . 192 1 [396[ 397]. Also macroscopic hydrodynamic-like con¬ 
tinuum equations for the macroscopic order parameters were obtained and analyzed 


Interestingly, when the above-mentioned density bands appear in one species, they can in¬ 
duce spatial heterogeneities in the other species via the inter-species coupling [395] . These 
results may be interesting for the dynamics of biofilm formation. Biofilms are surface- or 
interface-attached communities of microorganisms [ 398111399] , in nature usually composed of 
more than one species Eunmn]. Often at least part of the microorganisms in a him adopts 
a motile state ] 1TT(). 140211404j . Their collective behavior should be influenced by inter-species 
interactions. 

Other variants of the Vicsek model, such as polar particles of apolar alignment interac¬ 
tions [383 , 391.13931 1394 . 405 ] or nematic particles that randomly reverse their migration 
direction puniinH] were investigated. For the first kind of systems, density bands with par¬ 
ticle migration along the contour of the band were observed [3831.13911 [405] , in contrast to the 
above-mentioned density bands that migrate perpendicularly to their elongation direction. 


An example of apolar alignment interaction is depicted in Fig. 21 


A further characteristic feature of self-propelled particle systems has been outlined for the 
active nematic case [407] . We denote the average number of active particles in a subregion 
of the whole system by N and the fluctuations of this number by 8N. Typically a relation 
of the form ((5/V) 2 ) 1 ^ 2 oc is obtained, with £ = 0.5 in equilibrium systems. In the 
two-dimensional case of active nematics, an exponent £ = 1 was predicted [40711408] and 
later confirmed through simulations [ 406] . These large density fluctuations were termed 
“giant number fluctuations”. Likewise, in the case of active polar particles an elevated 
exponent £ > 0.5 was predicted [33B* 389] and confirmed in simulations [383-, 39~1 ,T(T5 . Such 
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large number fluctuations were also observed in experiments [36411366 1 #67]. As mentioned 
below, in systems featuring steric interactions, large number fluctuations can result from a 
cluster transition during which self-propelled particles mutually block their active migration 

pmmj . 

Finally, we note that two-dimensional systems of self-propelled particles can feature a 
long-ranged orientational order of their migration directions [ 388 . 1389] 1408] . This is in 
marked contrast to equilibrium systems [415j . Long-ranged order can emerge in the ac¬ 
tive case because information about the orientations is additionally propagated by the self- 
propelled motion of the particles. 


4-2. Self-propelled particles featuring steric interactions 

Originally, the point-like self-propelled particles in the Vicsek model pT7T] only interact 
via the alignment of their velocity orientations. In principle, this can lead to artifacts such 
as an unbounded growth of the local particle density. 

As a first step to solve this problem, modifications of the Vicsek approach were intro¬ 
duced. For example, the migration directions of particles that come too close to each other 
can be reversed to bound the density [ 416j . Or a pairwise alignment of the velocity vec¬ 
tors into opposite directions is induced for particles that are not well separated [417fl419j . 
Another route adjusts the migration direction according to an inter-particle potential [420]. 
If a finite intermediate distance is preferred, crystal-like arrangements can arise within the 
particle swarm [ 42'Oj . Such a situation can be interpreted as a combination of steric repulsive 
and cohesive attractive interactions EBB- 

In addition to steric interactions between the particles, also steric interactions with con¬ 
fining walls can be considered. An example is a channel geometry made by two parallel 
confining plates [418] , Using a functional form for the alignment interactions as in Eq. (15) 
in combination with discrete time steps, pattern formation was observed at intermediate 
overall particle densities in the channel [4T8] . On the one hand, lanes emerge along the 
channel direction that are either directly supported by the walls or “free-standing” in the 
interior [4T8) . see Fig. [22] 


(a). These lanes are “unidirectional” in the sense that all particles 
in the channel self-propel on average into the same direction. On the other hand, at slightly 
higher particle densities, active cluster crystals form that collectively migrate along the chan¬ 
nel BH. see Fig. i22\ (b). The mechanism behind the formation of these structures results 
from a combination of two ingredients: first discrete migration steps of the self-propelled 
particles, and second the possibility for overreactions in the velocity alignment. 

From an experimental point of view, dry systems of manifestly sterically interacting 
self-propelled particles can be realized by granular hoppers [364113671 l421H424j . For this 
purpose, the granular particles are typically set into motion by a vertical vibration of the 
horizontal substrate. Only the in-plane displacements are recorded. Both, the migrations 
and ordering of apolar [364. 121.14231 [424] and polar [365] rod-like objects were investigated, 
as well as the motion of polar disks [ 3661 367 . 425 ]. In the latter case, the steric interactions 
between the particles are isotropic due to the disk-like shape and do not explicitly provide 
an alignment interaction. Interestingly, the collective behavior of these vibrated polar disks 
could nevertheless remarkably well be mapped onto the Vicsek model |3S7j . 
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Figure 22: Pattern formation in a modified Vicsek model of self-propelled particles confined in a channel 
between two parallel horizontal plates [418] . The patterns appear at intermediate overall particle densities. 
As indicated by the arrows, the structures in both cases collectively migrate along the channel direction, 
(a) Increasing the overall particle density from below, first a state of unidirectional laning appears [418] : 
all particles on average migrate into the same direction, (b) At slightly higher particle densities, migrating 
cluster crystals emerge [ 4l8| . Each lattice point in this state is occupied by several self-propelled particles. 
The mechanism leading to the pattern formation is based on discrete migration steps of the self-propelled 
particles and the possibility of an overreaction in their velocity alignment. Channel walls are indicated by 
the gray bars and periodic boundary conditions are applied in the horizontal direction. 


A polar granular hopper features an anisotropic density distribution along its body to 
break the forward-backward symmetry. Under vibration it thus shows an imposed preferred 
migration direction [[365]. In contrast to that, for apolar hoppers the symmetry must be 
broken from outside by an additional horizontal mode of vibration [ 424] . or the symmetry 
must be broken spontaneously to perform steps of propagation. Vibrated apolar rods were 
observed to incline with respect to the substrate surface and to preferentially migrate into 
the inclination direction H2D 02 % The mechanism of spontaneous symmetry breaking 
was investigated in more detail for a dimer model particle [426] . In the case of vibrated 
polar disks, the polarity results from breaking the symmetry of contact with the vibrating 
substrate by attaching two asymmetric legs [367 ]. 

As revealed by these studies [4221,1426] , the “dry” friction between the particles and the 
substrate plays a major role for the propagation mechanism. The displacement statistics 
of “dry” objects on vibrated substrates were analyzed experimentally, also as a function of 
the substrate inclination [42714429] , Interestingly, the displacement distribution functions 
appear to have exponential tails [42714429] . Similar results were obtained for vibrated water 
droplets [4271 1430] . 

Typically the interactions with the substrate are modeled by a “dry” friction term of 
the Coulomb type [ 431] . Generally, a simple model equation for the one-dimensional motion 
of a single granular particle exposed to dry friction with the substrate and an additional 
dynamic (viscous) friction with the remaining environment can be written in the form of a 
Langevin equation as [43214437] 

dv v dx 

m— — —m - a(v)A + y(t), — = v. (16) 

dt t dt 

Here, m is the effective mass of the particle, v its (generally non-constant) velocity, and x its 
position. The first term on the right-hand side of the velocity equation denotes the dynamic 
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(viscous) friction with a relaxation time r. cr(v) returns the sign of v, i.e. a(v > 0) = +1, 
a(v = 0) = 0, and cr(v < 0) = —1. Thus the second term on the right-hand side includes 
the dry Coulomb friction of constant strength A. It implies that a moving particle is always 
slowed down by the same amount of deceleration, independently of its current speed. Finally, 
7 (i) describes a stochastic force, which may arise from fluctuations of the environment and 
stochastic noise in the vibration mechanism. Transferring Eq. (16) to the Fokker-Planck 


framework [T9TL11921 45? , 433 , T3Ti. 438] . one finds that the steady-state velocity distribution 
function features a cusp at v — 0 due to Coulomb friction [432H43611438, 439j . Interestingly, 
the irreversible part of the resulting Fokker-Planck equation can be mapped onto the case 
of a quantum mechanical harmonic oscillator with a pinning delta-potential at its center 
[436 ; 14404443] . Starting from a sharp initial state, the spatial particle distributions develop 
non-Gaussian tails on intermediate time-scales [436] . Multiple time-scales emerge in the 
evolution of the distribution functions due to Coulomb friction [ 436] . Nevertheless, despite 
such a non-Gaussian behavior, the mean-square displacement of the particles still grows 
linearly in time [1821 11831 436]. It will be an interesting problem for the future to study 
in more detail the effect of aspects of the friction mechanism on the collective behavior of 
granular hoppers. 

4-3. Active driving force instead of constant velocity 

In the previous sections, we reported about the success of the Vicsek model, which, 
despite its simplicity, could even reflect the collective behavior in an example system of 
granular hoppers [367] . One of the key assumptions in the Vicsek model is the constant self¬ 
propulsion speed of each individual particle [377 ]. This is an adequate simplification in a 
dilute system. However, it becomes problematic in dense systems of interacting particles. In 
particular, jammed situations where particles mutually block their ways cannot be described 
adequately in this model. 

It is then reasonable to switch to the picture of an active driving force [13611140911411114441 - 
454] instead of a constant self-propulsion velocity. In the simplest general case, the particles 
are now characterized by their position r, by their velocity v, and by a polar unit vector 
p that gives the current orientation of the active driving force. Here, situations can arise, 
in which the direction of the velocity v and the orientation p of the driving force are not 
parallel. The magnitude of the velocity v can change over time, whereas the strength of the 
active driving force FJ is usually kept constant for simplicity. 

In many cases, the motion of the particles is overdamped. This is true for instance 
for colloidal Janus particles |355H361j or bacteria [36811371] that self-propel in an aqueous 
environment at low Reynolds numbers [455] . Overdamped Brownian dynamics is appropriate 
to describe this kind of motion and the velocity variable becomes redundant. The structure 
of the rescaled dynamic equations for an isotropic self-propelled particle of constant strength 
Fd of its active driving force can be written as 


dt 


Fd P + £, 


dp , . 

a = (w+ v) x p 


(17) 


Here, the potential U contains for example the influence of external fields like gravity, the 
steric interactions with other particles, or the steric interactions with confining walls. The 
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angular velocity tjJ results from the torque on the particle for example due to external 
alignment by a magnetic field. Both, £ and 77, include the impact of stochastic fluctuations. 
So far, most of the corresponding studies have been restricted to two spatial dimensions. 

First, the stochastic motion of single isolated self-propelled particles was analyzed [ 356. 
1444] . In the long-time limit, the stochastic noise determines the statistical properties. Con¬ 
sequently, in the long-time regime, the motion appears diffusive. However, the magnitude of 
the corresponding overall diffusion constant can be significantly increased by self-propulsion 
when compared to the case of an analogous passive particle [ 356. 14441 I45'6| . In contrast 
to that, self-propulsion significantly alters an intermediate time regime. For example, 11011 - 
Gaussian displacement statistics can emerge at intermediate times in spite of a Gaussian 
form of the stochastic noise [444] . These results were confirmed by comparison to experimen¬ 
tal systems of self-propelling colloidal Janus particles |456j . Also the impact of an external 
flow field was studied [ 446 . 450] . 

Apart from that, Eqs. GD were generalized and investigated for the case of rod-like par¬ 
ticles, where the diffusion and friction matrices become uniaxially anisotropic |457 j. For gen¬ 
eral particle shapes, aligning torques can arise due to steric interactions between anisotropic 
particles or because of interactions with confining boundaries. These aligning torques enter 
the equations via the angular velocity cu. Furthermore, a torque that leads to a continuous 
reorientation of the polar direction p can arise from the mechanism of self-propulsion in 
combination with the particle shape or be actively generated. The stochastic motions of 
such circle swimmers were analyzed [ 4571, 458 j and compared with experiments [459], Self- 
organization in array-like vortex textures was found |460] . Moreover, it was demonstrated 


that in three spatial dimensions the asymmetry of biaxial self-propelled particles can lead 
to helical and even to superhelical trajectories [448] . 

Concerning the collective behavior of many interacting particles, we first focus on the 
case of spherical objects [ 36F, 40911411114451 447,1449] . Due to their isotropic shape, they 
do not feature a steric alignment interaction. This promotes the emergence of new truly 
non-equilibrium effects. A spherical particle that sclf-propels towards a hard wall is blocked 
in its motion and slowed down [461] , It can only escape again when the direction of active 
drive has reoriented away from the wall. The time scale for this process is set by rotational 
diffusion and thus indirectly by the temperature of the system. During this time, other 
particles can hit and additionally get blocked, if the density and self-propulsion speed are 
high enough. Through such a mutual blocking, see Fig. 23 (a), wall-supported clusters can 
emerg^J Interestingly, the analogous process is also observed in a self-supported way in the 
bulk [ 40911410] in two-dimensional systems, see Fig. 23 (b). Here, the blocking mechanism 
was analyzed in more detail [411] and also investigated experimentally [36 lj. A real phase 
separation into a clustered and a gas-like phase can be obtained as a steady state of the 
system [4491144311454] . For clarity, it is stressed that these clusters form in the absence of any 


3 This was demonstrated for rod-like particles [462] . In that case, after a certain time, the clusters can 
slip off the wall, supported by the reorienting steric interactions with the wall [4'6'1| . On the contrary, due to 
their isotropic steric interactions, self-propelled spheres are more effectively blocked by walls than rod-like 
particles [4611 . 
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Figure 23: Schematic illustration of the blocking mechanism of isotropically sterically interacting self- 
propelled particles. The self-driven particles are of spherical shape. Arrows indicate the direction of active 
drive which can reorient by rotational diffusion. A two-dimensional motion of the particles within the plane 
of the figure is assumed, (a) Two particles block each other, supported by a rigid wall, (b) Self-supported 
blocking can also occur in the bulk. If rotational diffusion is slow enough in comparison to the active drive 
and the density is high enough, both situations can serve as seeds for the formation of larger clusters. 

attractive forces and despite purely repulsive interactions between the particles. The clusters 
emerge solely from mutual steric blocking of the migration due to the non-equilibrium active 
drive. 

Continuum equations were derived that trace the phase separation into the gas-like and 
clustered phases back to a force imbalance arising from the mutual steric blocking [ 41 lj , and 
a weakly nonlinear stability analysis [4541] was carried out [ 465 j. In this way, the nature of 
the transition of phase separation could be investigated in more detail. Interestingly, the on¬ 
set of the clustering transition can be mapped onto the phase separation of passive particles 
described by the famous Cahn-Hilliard model [466 ]. In particular, this implies the existence 
of an effective free energy density that can characterize the non-equilibrium clustering tran¬ 
sition [465] . The analysis suggests a change of the clustering transition from continuous at 
high particle densities to discontinuous at lower particle densities. This change in the nature 
of the clustering transition was confirmed by numerical investigations [465j. Naturally, the 
clusters imply large fluctuations of the local particle density [409j. A reentrance of the fluid 
phase is observed as a function of the driving force [41 lj . 

At high densities, the active systems become jammed [ 445 , 449] or crystallize [ 4101, 447]. 
Active crystals were also observed experimentally mu, with attractive particle interactions 
typically being involved in the crystallization process [4671 168 . 

On the one hand, active crystals have been studied numerically by particle simulations 
! 110. 4l9l [4471. 451 . 1452 j . The self-organization of purely repulsively interacting actively 
driven particles into crystal-like structures was investigated [447]. Furthermore, the collec¬ 
tive behavior of self-propelled particles that are ordered in crystal-like arrangements and 
interact by harmonic-spring potentials has been analyzed in detail [4511 452] . In the latter 
case, there is no explicit alignment interaction between neighboring particles. Nevertheless, 
the elastic interactions between them channel the individual attempts of self-propulsion to¬ 
wards one global orientation of motion. Finally, the whole structure collectively migrates 
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into one common direction. 

On the other hand, a field-theoretic approach was introduced to study the collective be¬ 
havior of active crystals [ 469, 470] . It can be derived microscopically from classical dynamical 
density functional theory [47113479] applying appropriate assumptions. Two microscopic or¬ 
der parameter fields were used to characterize the state of the system composed of many 
self-driven particles. The first one is a particle-resolved density field -0 (r, t); the second one 
is a polar order parameter field P(r, t) that characterizes the locally preferred orientation of 
the active driving direction p, see Eqs. (17). Coupled dynamic equations describe the time 
evolution of these order parameter fields [ 469 . 470]. Resulting periodic structures are dis¬ 
played in Fig. [24} Concepts from two prominent continuum descriptions were unified in this 


approach. The first one is the phase field crystal model characterizing periodic modulations 
in the density field. This model was successfully used to reproduce solidification and crystal¬ 
lization phenomena in conventional passive crystalline structures on particle-resolved length 
and diffusive time scales [342113431134911473114761148011483] . The second prominent approach 
contained here stems from the macroscopic Toner-Tu model for the collective motion of 
non-crystallized self-propelled particles [38811389] , It can distinguish between systems that 
feature spontaneous directed self-propulsion, for example emerging from explicit alignment 
interactions, and systems without spontaneous alignment of their self-propulsion directions. 
Starting from initially disordered active systems, the field-theoretic approach shows that an 
active single crystal can form via the coarsening of multi-domain structures [469]. If a local 
alignment mechanism for the active driving directions of the particles exists, the crystal fi¬ 
nally collectively travels into one common migration direction. If there is no such alignment, 
the active crystal remains at rest for small active drive, see Fig. 24 (a). However, beyond a 
threshold value of the active drive, the crystal still finds a common migration direction and 
starts to collectively travel. This is similar to the above-mentioned results from the particle- 
resolved simulations of active crystalline structures [45111452] . With increasing active drive, a 
transition from traveling hexagonal to traveling rhombic, quadratic, and lamellar structures 
was observed [4~69] . Such a transition series is displayed in Fig. 24 (b)-(d). 

In addition to that, using the active phase field crystal theory [469] . a linear stability 
analysis of the traveling active single crystalline structures was performed [470]. This analy¬ 
sis is complicated by the fact that the density field in the crystalline ground state is already 
spatially modulated. It was found that the investigated collectively traveling active single 
crystals are linearly stable. The impact of hydrodynamic interactions on the stability of 
such active crystals was considered and is discussed in the next section. 

Turning now to non-spherical self-propelled particles, a well-studied example is given 
by self-propelled rods actively driven parallel to their long axis [ 4531 14601 14621 148413494] . 
The biological motivation arises for instance from experimental investigations on elongated 
bacteria crawling or gliding on a substrate [121EJ [3711484,495]. Another biological example 
are filaments actively driven by molecular motors. These motors can be fixed on a substrate 
at one end [37511496] or directly work between the filamental rods [37611497] . Non-biological 
artificial experimental systems can in particular be realized by vibrating polar granular rods 
[365]. In all these cases, the anisotropy of the steric interactions provides a mechanism of 
particle alignment. It turns out that the self-propulsion supports nematic ordering of rod- 
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Figure 24: Active crystals as modeled and displayed in Refs. 053 STD]. The particle-resolved density field 
ij)(r,t) is encoded by the color, where brighter color means higher density. Thin needles, pointing from 
the thicker to the thinner ends, indicate the local orientation and magnitude of the polar orientational 
order parameter field P(r, t). Thick arrows mark directions of collective motion. Panels from left to right 
were obtained for an increasing strength of the active drive. Here, explicit alignment mechanisms for the 
self-propulsion directions were not included. Therefore the hexagonal active crystal in panel (a) at low 
active drive remains at rest. Beyond a threshold magnitude of the active drive, however, this structure 
can nevertheless start to migrate, as shown for the traveling hexagonal crystal in panel (b). Increasing 
the active drive further, structural transitions to traveling rhombic and quadratic (c) as well as traveling 
lamellar textures (d) were observed. 


like particles by lowering the transition density from the isotropic to the nematic state 
This is in line with the prediction of long-ranged orientational order in two-dimensional self- 
propelled particle systems |388L 138911408| . Apart from that, various different collective states 
of migration have been revealed for interacting self-propelled rods [ 453. 462 . [4841.148811489. 
1491. 492], Among those are disordered, jammed, clustered, and swarming states, where in 
the latter case the rods organize in localized packets of common migration direction)^] Laning 
states were identified [ 4531. 489 . 491 . [492 ] that feature parallel stripes or lanes within which 
all particles propel into the same direction. However, neighboring lanes in these systems 
show antiparallel migration directions, see Fig. 25 for a schematic illustration. This is in 


contrast to the above-mentioned case of unidirectional laning [418] . Laning as in Fig. 25 


was 


previously observed in systems of externally driven particles, where for different particles a 
driving force into opposite directions was imposed from outside [498H505] . In the case of 
self-propelled particles, laning appears spontaneously without explicitly superimposing the 
opposite migration directions. Finally, even turbulent states were discovered in simulations, 
theory, and experiments on swimming bacteria [489] 1506] [507| , although these organisms 
swim in an environment of low Reynolds numbers. The source of this turbulence at low 
Reynolds numbers is the continuous energy input on the length scale of the individual self- 
driven particles. 

Apart from that, the behavior of active particles featuring various shapes more compli¬ 
cated than rod-like was studied [50811509] . In particular, forward-backward symmetries can 
be broken through the particle shape [508] . Identical shapes, but differently applied active 
driving forces can lead to phase separation [508] . Depending on the relative orientation of 


4 The notation in this context is not unambiguous. Sometimes polar swarms are also denoted as clusters, 
whereas sometimes non-polar clusters are denoted as swarms. 
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Figure 25: Schematic illustration of a laning state. Darker particles move upwards, whereas brighter 
particles migrate downwards. The particles organize themselves into neighboring lanes of opposite migration 
directions. Either the opposite migration directions of the individual particles are explicitly imposed from 
outside for particles driven by an external force; or the opposite migration directions emerge due to self¬ 
organization in systems of identical self-driven particles. In reality the lane thickness is not necessarily as 
regular as depicted here and can comprise different numbers of particles. 


the self-propulsion direction with respect to the particle shape, concave particles can self- 
organize into micro-rotors [508j. Finally, for binary mixtures of non-spherical active rotors, 
a phase separation into domains of clockwise and counter-clockwise rotations was identified 


4-4- Role of hydrodynamic interactions 

Many realizations of self-propelled particles exist in a liquid environment. Examples are 
artificial microswimmers as for instance the colloidal Janus particles addressed above [ 355k 
i361j , and swimming bacteria |368l - l37l] . Any self-propelling swimmer sets its surrounding 
fluid into motion. The motion of the other swimmers, which are suspended in the fluid, is 
naturally influenced by this fluid flow. In turn, also the motion of all the other swimmers acts 
back onto the first one. These fluid-mediated interactions are referred to as hydrodynamic 
interactions. They are long-ranged. 

To include hydrodynamic interactions, we must take into account the fluid flow. It is 
typically a very good approximation to consider the - often aqueous - surrounding liquid as 
incompressible, i.e. to set its mass density p equal to a constant. Then the central equation 
of motion in fluid dynamics that determines the fluid flow, the Navier-Stokes equation [ 510i| . 
can be written in the form 

P + P [ v ( r : t) ■ V] v(r, t) = — Vp(r, t) + ?/V 2 v(r, t) + f(r, t). (18) 


It describes the influence of volume force densities acting on a fluid element at position r 
and time t. v(r,t) is the resulting flow velocity field of the fluid at position r and time t. 
The left-hand side of Eq. (18) gives the inertial contribution; here the second term expresses 

44 


















the fact that the actual motion of the fluid elements must be taken into account to correctly 
determine the change of momentum at a spatial position r. On the right-hand side, the first 
term denotes the force density due to gradients in the pressure held p(r,t), the second term 
includes dissipation with rj the viscosity of the fluid, and the last term takes into account 
any additional force density held f(r, t) acting on the fluid. In our case, f(r, i) contains the 
force densities that each swimmer exerts on its liquid environment. 

Rescaling all lengths by a characteristic length scale L , all velocities by a characteristic 
velocity V, time by L/V, pressure by pV 2 , and force density by pV 2 /L } Eq. (18) takes the 
form 

— + [v 7 ■ V 7 ] V 7 = -vy + Re-W + f 7 , (19) 

with space and time dependencies not marked explicitly any more. We find 


Re = 


pLV 

V 


( 20 ) 


for the famous dimensionless Reynolds number. Typical microswimmers in the form of 
colloidal Janus particles and swimming bacteria have dimensions L in the range between 
10 nm and 10 pin. A characteristic swimming speed would for example be 10 pm/s. Thus, 


we can see from Eq. (20) that their swimming in an aqueous environment usually occurs at 
low Reynolds numbers Re < 1 [355] . 


We can infer the consequences of this estimate from Eq. (19). For low Reynolds numbers, 


the inertial terms on the left-hand side can be neglected when compared to the dissipative 
contribution on the right-hand side. Hydrodynamics at low Reynolds numbers is thus de¬ 
scribed by Stokes’ equation [ISTTi]: 


0 = -V 7 p 7 + Re-W + f 7 . 


( 21 ) 


This equation leads to some peculiarities for swimming at low Reynolds numbers that are 
markedly different from our every-day experiences. In particular, the equation is linear 
in the velocity held and instant in time. Consequently, for balanced pressure gradients, 
how immediately stops when no further force density f 7 is applied to the fluid. This has 
qualitative consequences when we compare to the situation at higher Reynolds numbers. 

At high enough Reynolds numbers, a swimmer can move by reciprocal shape changes, 
see also Fig. 26: hrst a quick shape change is performed in the form of a power stroke; then 


a slow recovery stroke follows in the exactly inverse way of the power stroke, only that it is 
performed significantly more slowly. So the central difference between the two strokes lies in 


the speed of execution. In this situation, the inertial terms on the left-hand side of Eq. (18) 


can lead to a symmetry breaking and net motion. However, reciprocal shape changes of the 
swimmer cannot lead to such net motion at low Reynolds numbers. The symmetry-breaking 


inertial terms are missing in the Stokes equation Eq. (21). These facts are often referred to 
as the “scallop theorem” [455]. 

Furthermore, in the absence of pressure gradients, the swimmer motion completely stops 
together with the fluid flow v 7 as soon as the swimmer does not apply any force density 
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Figure 26: Illustration of the qualitatively different situation for swimming at high and low Reynolds number 
(“scallop theorem”) I I'm! . The simple swimmer (black) performs a reciprocal shape change, i.e. the shape 
change during the power stroke is just the inverse of the shape change during the recovery stroke. Only the 
speed of shape change is significantly higher during the power stroke than during the recovery stroke. At 
high Reynolds numbers this leads to a net displacement after the cycle is completed due to inertial effects. 
At low Reynolds number, inertial effects are negligible and a net displacement cannot be achieved in this 
way. 


f' to the fluid any longer. Apart from that, the sum of all instant forces exerted by a low 
Reynolds number swimmer on the surrounding fluid vanishes because each stroke that it 
performs is balanced by a counter-acting drag force due to the motion of its body. It can 
be shown that, with increasing distance from such a “force-free” swimmer, the induced flow 
field decays significantly more quickly than for an object that exerts a net force on the 
surrounding fluid E2|. Furthermore, the swimmer must provide a suitable mechanism to 
break the symmetry and achieve a net forward motion [ 513j . 

Different routes to reach this goal were pointed out. One example are non-reciprocal 
cycles of contraction and expansion of a straight three-linked-sphere swimmer [5141 515J as 
displayed in Fig. 27 (a). Another example are volume changes of the spheres of a two-linked- 
sphere swimmer during contraction and expansion [516] . see Fig. |27| (b). Apart from that, 
the rigid structure of a swimming object can already break the symmetry by itself. Ideal 
candidates for this purpose are helical objects that have a certain handedness [ 369 : l5T8h 
1523] . And as mentioned above, Janus particles intrinsically break the forward-backward 
symmetry, which can be exploited to start a net-propulsion mechanism [3551 - 136ll 1363 ]. 

One swimming strategy for a swimmer at low Reynolds numbers is to exploit the drag 
resulting from its self-induced fluid flow. To understand this principle, it is most illustrative 
to consider two force centers as shown in Fig. 27 (c). In the upper case, the forces applied to 


the fluid are oriented in a way that the swimmer pushes the fluid outwards. It is therefore 
called a pusher. The bottom case shows the opposite situation, in which the swimmer pulls 
the fluid inwards. Thus the swimmer is referred to as a puller. If the swimmer body is 
located asymmetrically with respect to the centers of force acting on the surrounding fluid, 
it will experience a drag from the self-induced fluid flow [ 5171 524 , 525 ]. Thus the swimmer 
will show a net propulsion. This propulsion principle can be generalized to more realistic 
swimmer shapes. 
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Figure 27: Illustration of different self-propulsion mechanisms for swimming at low Reynolds numbers. Panel 
(a) shows the three-linked-sphere swimmer suggested in Ref. |514j . Its deformation cycle (snapshots from 
top to bottom) is nonreciprocal. Due to hydrodynamic interactions, contraction or expansion of one bond 
leads to a different net displacement depending on whether the other bond is in the contracted or expanded 
state. The vertical bar is added for clarity to highlight that a net propulsion to the right has occurred during 
the cycle. Panel (b) depicts the nonreciprocal deformation cycle (again snapshots from top to bottom) of 
the two-linked-sphere swimmer introduced in Ref. [ 15161 . In the completely contracted or expanded state 
of the bond, the swimmer exchanges the volume of the two spheres, for example by pumping the contents 
from one side to the other through the bond. Since the viscous drag is higher for the sphere of larger radius, 
the swimmer features a net propulsion to the right (again the vertical bar is added to make the net motion 
visible). Finally, panel (c) displays the minimal model swimmer suggested in Ref. 15171 . The outer spheres 
are small and experience a negligible viscous drag by the fluid. However, they represent centers of forces 
acting on the fluid and setting it into motion (indicated by the bright arrows). Due to the antiparallel 
alignment of the forces, the total applied force vanishes. In the upper case, the swimmer pushes the fluid 
outwards, which is why it is referred to as a “pusher”, while in the lower case, it pulls the fluid inwards and 
would be called a “puller”. The main swimmer body is asymmetrically located with respect to the force 
centers. Thus it experiences a net viscous drag by the self-induced fluid flow. 


We now turn to the situation of more than one swimmer suspended in the fluid. Then 
one of them feels the drag due to the resulting fluid flow induced by all other swimmers, 
and vice versa. The impact that these hydrodynamic interactions between active swimmers 
has on their net interactions and on their collective behavior is still under investigation. 
For the three-linked-sphere swimmer mentioned above [514j . the pairwise hydrodynamic 
interaction was analyzed explicitly |52b| . A complicated dependence on the relative distance, 
orientation, and phase of the shape change was obtained. Another study focused on the 
averaged flow fields around circle swimmers [527] . For particles that propel by pushing 
the fluid, a repulsive interaction was found on average, whereas for swimmers that pull the 
fluid, the average interaction was attractive. A pair of counterrotating pullers performs a net 
translational motion [ 527( i 528] . Straight-swimming pullers that were forced to swim parallel 
to each other were observed to repel each other [ 52911530] . whereas pushers in the same set¬ 
up were found to attract each other [506 , 529) . However, when released from the constraint 
of parallel swimming, no permanently bound state due to hydrodynamic interaction was 
found in spite of the initially attractive configuration [529} . 

Several cases of motion in a liquid environment were identified in which hydrodynamic 
effects seem not to play a crucial role. For bacterial cells, it was demonstrated that generally 
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cell-cell interactions and ccll-snrface interactions are dominated by orientational diffusion, 
steric interactions, and lubrication forces, and that hydrodynamic interactions play a mi¬ 
nor role put Apart from that, different collective effects observed experimentally for mi¬ 
croswimmers were successfully reproduced in simulations without including hydrodynamic 
interactions. Examples are the previously addressed observations of clusters and turbulent 

states mm- 

Nevertheless, there are many situations in which hydrodynamic interactions were shown 
to be important, both for the single-swimmer and for the collective behavior. For instance, 
hydrodynamic interactions with confining walls can modify the migration properties [5311 - 
l536j . As demonstrated numerically, the nature of the propulsion mechanism can determine 
whether active microswimmers preferentially orient towards a confining surface, or whether 
they turn away and leave the wall j[535] . Under narrow confinement between two walls, 
the propulsion mechanism can support orientations towards one of the walls or orientations 
parallel to their surfaces, which supports blocking or planar swimming motion, respectively 


[361 , 53511536] . Apart from that, a speed-up of microswimmers through hydrodynamic in¬ 
teractions with a flexible confining tube was predicted ESI- Microswimmers conhned in 
a harmonic trap were found to orientationally order due to hydrodynamic interactions and 
induce a net fluid flow [538]. Moreover, in experiments, a transient capturing of microswim¬ 
mers in circular trajectories around passive colloidal spheres was observed [539J . 

Generally, hydrodynamic interactions can induce the formation of vortices and swirls 
[Sffij- In that sense, hydrodynamic interactions can destroy long-ranged orientationally or¬ 
dered collective motion, but also more localized orientational order in swarms [54lj . Likewise, 
in the case of two-dimensional active traveling crystals, destabilizing effects of hydrodynamic 
interactions were observed numerically [470]. The self-propelling particles forming the active 
crystal were considered to migrate on a substrate covered by a thin fluid him. As a con¬ 
sequence of the resulting how helds, the traveling single crystal can break up into different 
domains, or the crystalline lattice structure can vanish altogether. The stability of different 
lattice structures under hydrodynamic interactions was analyzed [542] . 

Recently, it became clear from theoretical and numerical investigations that hydrody¬ 
namic far-held and near-held interactions can differ in their influence on the orientational 
order and stability of the suspension. Naturally, far-held interactions prevail in dilute sus¬ 
pensions whereas near-held interactions become important in concentrated solutions of mi¬ 
croswimmers. Most of these studies were performed using yet another swimmer model called 
“squirmer” [543H545] : a non-vanishing how held is prescribed as a boundary condition for 
the surrounding fluid on the surface of each swimmer; the swimmers themselves are often 
considered as rigid objects, but the non-vanishing surface how drives them through the liquid 
environment. 

Following these lines, we hrst consider dilute suspensions of swimmers that solely interact 
hydrodynamically. For those, it was demonstrated that orientationally ordered states are al¬ 
ways unstable [546R548] . Furthermore, for elongated swimmers featuring a puller propulsion 
mechanism, it was found that isotropic and spatially homogeneous suspensions are stable 
|546fi548j . In contrast to that, homogeneous isotropic suspensions of elongated pushers are 

48 


[ 535 . 1536 ], Naturally, this influences the formation of the clusters mentioned in Sec. 4.3 







































unstable and show an interesting nonlinear dynamics including enhanced mixing [54611548] . 

When concentrated suspensions or the close proximity of confining boundaries are ad¬ 
dressed, hydrodynamic near-field interactions become important, ft is then often mandatory 
to explicitly resolve the induced flow fields in the vicinity of the swimmers. For this purpose, 
squirmer models are ideal candidates 1521) . 15351153611549. 550] . In the case of concentrated 
suspensions of spherical squirmers in the bulk, net polar orientational ordering of the swim¬ 
ming directions due to hydrodynamic interactions was numerically observed [ 54911550 ] . This 
is in contrast to the above-mentioned results for dilute swimmer suspensions [54611548] . For 
pullers the degree of order was higher than for pushers [ 549 . 550 ], Likewise, as indicated 
above, hydrodynamic near-field interactions between two confining walls determine the rel¬ 
ative orientation with respect to the boundaries [ 535 ; , 1536] . An enhanced cluster formation 
results for pusher squirmers when compared to pullers 535 [ 536 ] . 

Finally, hydrodynamic interactions can lead to synchronization. This particularly applies 
when active rotors or driven filaments are considered [55111567] . It was summarized that rotor 
pairs can only synchronize if the system is not symmetric under exchanging the two rotors 
[561!, 568] . Carpets of driven hydrodynamically coordinated active rotors were suggested 
to be used to pump or mix fluid in microfluidic devices [557]. Furthermore, carpets of 
actively driven deformable filaments are found on the surfaces of certain bacteria and can be 
used for self-propulsion 1171. 553 . 562] . Their coordinated motion can provide a swimming 
mechanism for these cells at low Reynolds numbers. 

4-5. Deformable self-propelled particles 

On our way of increasing complexity in active particle systems we now add a further 
degree of complication. Several self-propelled objects are not rigid. They deform as it is 
observed for instance for certain cells that crawl on substrates [372, 569. 1570 ], Particularly 
obvious examples are self-propelled droplets on interfaces [57111573] or in bulk fluid [57411576] . 
These droplets maintain chemical reactions [57411576] that lead to concentration gradients 
along their surfaces. A similar mechanism works by depositing chemicals from their inside to 
the outside [57114573] . In effect, all these processes induce and maintain gradients of surface 
tension along their surfaces [ 5731 15 7514578] . These lead to stress gradients along the surfaces. 
On rigid substrates, such stress gradients on the droplet surface can directly push or pull 
on the droplet and lead to a net drive. In fluid environments, the surface gradients can 
induce convective flows on the inside and outside of the droplet surface. These convective 
flows in turn propel the droplets [575 , 15771 1578] . If no special action is taken from the 
outside, the droplet motion has to start by spontaneous symmetry breaking; for instance 
by a sufficiently strong concentration fluctuation along its surface. In the ideal case, the 
concentration fluctuation leading to the net motion is maintained in a self-supported way by 
the resulting motion [ 579 ] . for example when the concentration gradients are enhanced by 
the induced fluid flows. It was demonstrated that self-propelled droplets can be relatively 
robust deformable objects. In an experiment, they were observed to deform and squeeze 
themselves through steric barriers [ 573 ] and recover afterwards. 

To study the impact of deformability on the particle motion, shape deviations from a 
spherical undeformed state can be taken into account in a systematic way [ 580] . The lowest 
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Figure 28: A self-propelled particle that deviates from its spherical ground state by axially symmetric 
elliptic deformations. The elliptically deformed state is described by the degree of deformation s and by the 
orientation of the resulting ellipsoid. Since there is no deformational forward-backward asymmetry along 
the symmetry axis of the ellipsoid, we indicate its orientation by a double-headed arrow referred to as n. 
The self-propulsion contributes to the velocity v of the particle. Although representing separate degrees of 
freedom for each particle, the velocity v and the deformation tensor S = s(nn — I/d) are generally coupled 
to each other, see Eqs. (22) and (23). 


order deformations are of axially symmetric elliptic shape. Such elliptic deformations can be 
described by a symmetric traceless tensor S of second rank [ 581 ]. It is formally equivalent to 
the order parameter tensor introduced in Eq. (J2| to characterize nematic liquid crystalline 
phases [5j. Now, however, s parameterizes the degree of deformation, while n represents 
the orientation of the symmetry axis of the elliptic shape changes, see also Fig. 28. The 


relation to the nematic order parameter tensor comes from the fact that axially symmetric 


elliptic deformations are forward-backward symmetric. In Fig. 28 this is indicated by the 
double-headed arrow that marks the symmetry axis of the ellipsoid. 

Coupled dynamic equations between the elliptic deformation tensor S and the velocity 
vector v of the particle were derived on symmetry grounds and evaluated [ 582] . To lowest 
order in the coupling terms, they read 


( 22 ) 

(23) 


dv 2 a 

—— = yv — v v — dS'v, 

at 

dS Q , i f 1 2t\ 

— = —kS + 6 vv - v I). 

at \ a 


In passive systems, the coefficient 7 would be negative and the term would correspond 
to linear viscous friction with the environment. On the contrary, in active systems, this 
coefficient can become positive, 7 > 0. Then, the first two terms on the right-hand side of 


Eq. (22) can lead to a stationary solution of non-vanishing velocity v/0, i.e. self-propulsion. 


An analogous approach can be found in the famous macroscopic continuum characterization 
of flocks of self-propelled particles by Toner and Tu [ 388 . 1389] . The orientation of the 
resulting non-vanishing velocity vector v is not fixed a priori and results from spontaneous 


breaking of the continuous rotational symmetry. In the second equation, i.e. in Eq. (23), the 
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coefficient n > 0 expresses that the shape tends to relax back to the nndeformed ground state 
for vanishing velocity v. As before, vv represents a dyadic product, d the dimensionality 
of the system, and I the unity matrix. Most importantly, the terms with the coefficients 
a and b include the leading-order coupling terms between the dynamic variables v and S. 
These contributions are to a big extent responsible for the peculiar behavior of deformable 
self-propelled particles as it is briefly summarized in the following. 

In two spatial dimensions, the deformability can induce a bifurcation from straight to 
circular motion with increasing active drive [582] . Also quasiperiodic motions as well as 
array-like trajectories featuring rectangular edges were observed [ 580] . Applying an exter¬ 
nal electric held to polarizable particles or an external gravitational held, further types of 
trajectories such as different kinds of zig-zag motions and cycloidal motions appear [ 583] . In 
the absence of external helds but in three spatial dimensions, additional helical trajectories 
are obtained [584 ; 585 ]. 

The collective behavior in two spatial dimensions was studied extensively [582; 586]. 
First, a global coupling was introduced that forces the deformed particles to align along their 
globally averaged deformation axis [ 582] , This coupling can lead to chaotic motion [ 58211586] . 
After that, instead of the global coupling, a pairwise alignment interaction between deformed 
particles was imposed [ 587 . 588] . In addition to that, as a steric interaction, a soft pairwise 
repulsive Gaussian potential was applied. Starting from random initial conditions at low 
or vanishing orientational noise amplitudes, the particles collectively order their migration 
directions as a function of time as expected [ 377 ] . All particles then collectively propel into 
the same direction. Moreover, they also tend to locally positionally order in a hexagonal 
way. Interestingly, it was observed that under compression of the system, implying an 
increase in the particle density, the ordered collective migration breaks down at a critical 
density [587 l 1588] . The system becomes disordered and fluid-like. This is just the opposite 
effect as to that expected from the Vicsek model, in which the orientational order in the 
motion increases continuously with the particle density [ 577 ] . In a later study that used 
an anisotropic pairwise repulsive Gaussian potential, it was concluded that the disorder 
transition from a collectively traveling hexagonal crystal to a fluid-like state is enabled by 
the soft Gaussian interaction potential [589]. The effect is well known in corresponding 
equilibrium systems. At low densities, these equilibrium systems are in a fluid phase, they 
crystallize at intermediate densities, but they return to a fluid state at high densities [59(Jh 
A central ingredient for this reentrant fluidization is the boundedness of the Gaussian 


interaction potential, even when two particles are placed at the same position [590] . It is 
concluded that the analogous effect is observed in active systems of self-propelled particles 
interacting sterically by a pairwise Gaussian potential [ 589] and can be supported by their 
deformability [587H589J . 

Apart from that, mutual deformations of the particles due to steric interactions between 
them were considered [589]. This process can provide an implicit alignment mechanism for 
their self-propulsion directions as illustrated in Fig. 29 for two particles. In many-particle 
systems, hexagonal lattices of self-propelled particles featuring a common migration direc¬ 
tion can emerge that represent active collectively traveling crystals. Introducing a cut-off for 
the steric Gaussian interaction potential, traveling crystals of rectangular lattice structure 
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Figure 29 : Aligning steric interactions between deformable self-propelled particles [555]. Depicted is a time 
series from bottom to top with dotted lines indicating the particle paths. Two deformable self-propelled 
particles migrate towards a common point. At this point, their active drive pushes them together. Due to 
steric interactions, this leads to deformations. Eqs. ( 22 ) and ( 23 ) imply that the deformations can reorient 
the self-propulsion directions, here towards the common axis of elongation of the particles. The outgoing 
angle between the particles is smaller than the ingoing angle, which means that an effective alignment 
mechanism is at work. 


are obtained. Likewise, active resting cluster crystals with more than one particle on each 
lattice site are observed, as they were previously reported for equilibrium systems |594[ 1595] . 
In cases of strong particle deformations, laning patterns can appear from randomly initial¬ 
ized systems [589] . Adding orientational noise, the Vicsek transition [377] from collectively 
ordered to disordered motion is recovered with increasing noise amplitude [589] . 

Finally, the formation of propagating fronts was reported for two different cases. First, 
large systems of repulsive Gaussian interactions and reentrant fluidity become inhomoge¬ 
neous at the reentrance density [38711588 ]. They show regions of collectively ordered and 
regions of disordered motion. In contrast to the traveling density bands observed in the 
Vicsek model [3791 1380 . 3831 1395] , here a front of disordered motion propagates into the 
region of collectively ordered migration. Furthermore, the disordered region has a higher 
density than the area of ordered collective motion. Second, and in analogy to the results 
from the Vicsek model [379, 380, 383 . 3951 . traveling high-density bands of ordered collective 
motion emerge in systems of density-dependent active drive [386] . The active drive in this 
case is set to increase with density. Such traveling density bands were shown to survive and 
penetrate through each other under repeated head-on collisions [386] . 

Apart from that, the model was complexified in two different directions. On the one 
hand, the next-higher mode of deformation was included [ 580] . It can be described by a 
third-rank tensor that is familiar from the study of liquid crystals of tetrahedratic 

order [ 59611597] . This uneven mode of deformation breaks the forward-backward symmetry 


52 






































Figure 30: Elongation of a deformable particle by an externally imposed shear flow [ 1603 ]. Initially spherical 
deformable particles like droplets (left) are deformed due to the elongational component of the flow (right). 
In addition to that, the rotational component of the shear flow tends to induce further dynamical processes, 


see Figs. 11 and 17 


of the particle shape. In this case, a zig-zag trajectory and chaotic motions can be obtained 
already for a single particle in the absence of an external field [580], The description was 
further extended to include the lowest four modes of deformation |598j . On the other hand, 


in addition to self-propulsion and deformability, an active rotational (spinning) motion was 
experimentally observed [ 599j and theoretically considered [60 0 ]. Different kinds of circular 
and quasi-periodic orbital trajectories result from this modification [600] . Furthermore, 
period-doubling and chaotic behavior is observed in the trajectories and in the state variables 
such as in the degree of deformation [50D] . Helical types of motion are found in three spatial 
dimensions, among them also one with a superhelical trajectory [ 60 lj . 

In a next step, the two-dimensional dynamics of a deformable self-propelled particle in 
a swirl flow [602] and in a planar linear shear flow were analyzed [6j03J. When exposed to 
shear, the elongational part of the flow tends to additionally deform the particle as depicted 
in Fig. [30] Moreover, in the case of the planar linear shear flow, the rotational contribution 


of the flow field tends to reorient the particle axes and can complexify the dynamic behavior. 
For this situation, an active straight motion, a winding motion, as well as different types of 
cycloidal motion were predicted [ 603] . An interesting issue for future studies is to find out, 
whether these kinds of motion represent the active analogues to the flow-aligned or tank- 


treading, wagging or swinging, and tumbling modes, respectively, displayed in Figs. 11 and 


17 


Experimentally, the shear geometry could be realized by placing self-propelling droplets 
on the surface of a liquid [57211604] that is confined between two parallel counterpropagating 

walls. 

The situation gets even more complex, when an additional active rotational spinning of 
the particle is considered in the shear flow [ 603] . In this case, active spinning with the same 
sense of rotation as the externally imposed flow field must be distinguished from active 
spinning with the opposite sense. When the active and the externally imposed rotations 
counteract and balance each other, a kind of active straight trajectory is recovered. Other 
types of observed dynamics are periodic, quasi-periodic, different kinds of cycloidal and 
undulated cycloidal, different kinds of winding, as well as different kinds of chaotic motions. 

To summarize the situation, we can say that the additional degrees of freedom arising 
from the deformability lead to quite complex dynamics, even for isolated active particles. 
Combining deformability, active propulsion, active rotation, and the various types of external 
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influences - in our example the imposed shear flow - opens the way to a nearly unbounded 
number of dynamical states. 

4-6. Collective behavior of animals 

At the end of this review, we now reach the uppermost spring of our ladder of complexity 
in the field of active systems. By construction this is a highly non-trivial topic and we can 
only touch the surface of it in the present framework. The main question that needs to be 
answered in this context is how much the characterization of each of the following systems can 
be simplified. The physicist’s goal is naturally to map the animal behavior onto a minimum 
model that captures the central ingredients but can still be evaluated efficiently. How will 
this work for living creatures that have their own mind and make their own decisions? 

Already for microorganisms the situation is quite difficult [353]. For example, reversals 
or changes in their migration directions are often observed [ 3681 5551160514608] . The latter 
can lead to a type of run-and-tumble motion [ 368 . 1463115551160714609] . On the one hand, 
this behavior may result from the complicated propulsion mechanism [609] . On the other 
hand, it can be the manifestation of an “active decision” iin!. If the change in migration 
direction occurs in response to an external stimulus, the reaction is generally referred to as 
“taxis” or “tactic” behavior. For instance, microorganisms can react in a chemotactic way 
when they detect chemicals that indicate nutrition, or they use chemicals for communication 
between each other without physical contact [ 368 . 4001 16071, 161114614] . Further examples 
concern the reaction to light (phototaxis) [370] 161511616] . gravitation (gravitaxis) [6171 618]. 
flow fields (rheotaxis) [619] . adhesion gradients (haptotaxis) [62011621] . and other external 
stimuli. Nevertheless, aspects of the collective behavior can often be described by simple 
particle descriptions in the form of extended Vicsek approaches [ 373 , 1416] or by minimum 
continuum models [489114951150711622] . 

Grasshoppers (or locusts) form a class of insects the collective dynamics of which has 
attracted attention for at least one reason: cannibalism is observed for these insects and 
seems to influence their individual and collective behavior [62314625] . On the one hand, an 
individual tries to rush after grasshoppers ahead of itself. On the other hand, it tries to flee 
from individuals behind itself. This was demonstrated by restricting the visual field through 
paintings on part of their eyes [626 ]. A restricted ability to detect approaching members in 
a group reduces motion and promotes cannibalism [626 ], which indicates that cannibalism 
supports migration. The situation was modeled in a pursuit-escape approach by finite-size 
random walkers that experience friction with the substrate and social interactions [ 627, '628]: 
other individuals moving away in the front are hunted after, whereas they are fled from if 
they approach from behind. In particular, the pursuit interaction tends to increase the 
orientational order in the collective motion. It increases the mean migration speed, which 
is also the case for the escape interaction at higher densities. A further study analyzed the 
radial distribution of other grasshoppers around one individual [629] . It revealed an isotropic 
shape in contrast to the anisotropic one obtained from minimum pursuit-escape approaches. 
The situation might still be more complex. 

Another group of animals repeatedly studied in the context of collective motion are 
fish. Data can be obtained by video-tracking the individuals of a swarm in artificial water 
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tanks [6301 - 1636] . It was observed that fish swarms were of smaller size when food sources 
were detected [637] and of larger size in reaction to predator attacks [637] . The swarms 
are generally elongated and denser at the front than at the back [6341I638l| . Increasing the 
mean density of fish in a container, a transition to cooperative collective motion could be 
detected at a threshold density [631 , 1639] similarly to the predictions of the Vicsek model 
[377] . However, the interaction rules deduced for individuals appear to be quite different 
from those assumed in the Vicsek model: direct orientational alignment seems to be weak 
[ 633, 634 ]. interactions are restricted to the one single nearest neighbor [633] or at least are 
not a simple average over all pairwise interactions [634], and variations of swimming speed 
are important [ 633 , 634 ], In a recent study, three modes of collective swimming motion 
were identified [636]: a swarming state of relatively low order, a polarized state, and a 
milling state of high global rotation. Transitions were frequently observed. These forms of 
collective behavior should be compared to the motion of single isolated fish. In the latter case, 
relatively continuous trajectories without sharp bends or kinks were observed, corresponding 
to persistently turning walkers rather than conventional random walks [63211635] . 

While the swimming motion of small fish can still be effectively confined to a quasi 
two-dimensional geometry in the lab, the study of fish and bird swarms in nature must 
typically consider three-dimensional textures. An exception are linear string-like and V- 
shaped arrangements (skeins) of, e.g., migrating geese, the structure and dynamics of which 
was recorded and analyzed [ 6401 1641 ]. To empirically study real three-dimensional flocks 
of birds, however, new tools had to be developed to effectively reconstruct the individual 
positions within the swarm from recorded two-dimensional image data [64211644] . In contrast 
to the observations made for fish, starling swarms were denser at their outer boundary than 
in the center [ 644 , 1645] without a systematic front-back asymmetry [ 645j . A correlation 
between the density and number of individuals in the swarm could not be identified [ 645] . 
On average, the probability to find a nearest-neighboring bird in the direction of motion was 
markedly depleted compared to other directions [645 , 1646] . Maybe most importantly, the 
interactions between the individuals were observed to follow topological rather than metric 
rules [64611647] : the interactions occur with a fixed number of nearest neighbors instead of all 
neighbors within a fixed interaction distance. Qualitative differences arise from this subtlety, 
for example the swarms are more cohesive under predator attacks [646] . Furthermore, in the 
famous Vicsek model [47?] . a change to non-metric topological interactions renders the order- 
disorder transition for collective motion continuous, without the emergence of density bands 
[392], This example demonstrates the importance of intensified recording and comparing 
to empirical data - a task that is highly non-trivial for such complex systems as animal 
swarms. 
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5. Conclusions 

In the above, we took a tour through various subtopics concerning the physics of soft 
matter out of its equilibrium ground state. The held of soft matter has grown far too broad 
for all its aspects to be touched in such a brief account. Only selected issues could be 
included. We were guided by two central aspects characteristic for soft matter systems: 
they typically show large responses to external stimuli; and they are easily driven out of 
equilibrium. The degree of non-equilibrium served to order the different topics, where in 
each case we considered systems of increasing complexity. 

Static external fields can switch the state of a system from its initial equilibrium ground 
state to a new static equilibrium state. We considered external electric fields that reorient the 
director in cells of low-molecular-weight liquid crystals as well as in swollen or prestretched 
liquid crystalline elastomers. In the latter materials, also mechanical holds can be applied 
to change the director orientation. Furthermore, the electric helds can induce mechanical 
deformations in swollen liquid crystalline elastomers, which makes these materials candidates 
for the use as soft actuators. The same is true for ferrogels or magnetic elastomers when 
they are exposed to external magnetic helds. Apart from that, the external helds can be 
used to tune the mechanical properties. For example, the elastic moduli can be reversibly 
adjusted in a non-invasive way by applying an external held. 

After that, we addressed externally imposed shear hows that lead to steady and dynamic 
states of motion. In such shear hows, the director of nematic liquid crystals can be observed 
to “how align”, i.e. to take a steady inclination angle with respect to the how within the shear 
plane; it can oscillate in a “wagging” motion; or it can describe continuous full “tumbling” 
rotations. In certain cases, these states can also be observed successively with decreasing 
shear rate. Orientational effects under shear are further found in bulk-filling, but spatially 
periodically modulated systems. We considered the example of micro-phase-separated block 
copolymer melts or solutions. Especially under large-amplitude oscillatory shear, reorienta¬ 
tions of the structures with respect to the shear directions occur as a function of the shear 
rate and amplitude. Furthermore, for localized objects such as vesicles that are composed 
of closed bilayer membranes a picture similar to the one for the director of nematic liquid 
crystals emerges under steady shear how. The steady state of constant inclination angle is 
called “tank-treading”, with the bilayer membrane running around the interior of the vesi¬ 
cle; an oscillating inclination angle is referred to as “swinging”, “trembling”, or “vacillating 
breathing”; and, again, full rotations occur in the “tumbling” mode. A non-vanishing vis¬ 
cosity contrast between the liquid on the inside and on the outside of the vesicle is necessary 
to obtain these states successively with decreasing shear rate. 

Finally, we turned to systems composed of constituents that themselves are considered to 
be active. They feature a mechanism of self-propulsion and possibly also of active rotations. 
Our focus was on the collective non-equilibrium states arising from the interactions between 
these constituents. Examples are ordered collective motion, traveling density bands, laning, 
resting and traveling crystals, or non-equilibrium clustering. We increased the complexity 
from point-like particles to finitely-sized objects including steric interactions; from particles 
propelling with constant velocity to those featuring a constant strength of active drive; 
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from dry systems to those interacting hydrodynamically; from rigid particles to those being 
deformable; and, eventually, by briefly addressing recent studies on living creatures. Quite 
remarkably, basic aspects, such as the formation of flocks of individuals that order their 
migration directions at high enough density, can be found on all levels of complexity. This 
supports the idea that it is possible to develop basic underlying concepts and systematic 
methods to describe also such kinds of genuinely non-equilibrium systems. 

As a concluding remark, we would like to append a central argument in favor of this held. 
Although often emerging from topics of daily life or biology, many of its subareas are rela¬ 
tively young. The boundaries are rapidly being pushed forward, thus continuously providing 
new interesting problems of research. Far from equilibrium, i.e. inherently for active systems, 
the tools of understanding still need to be established, although often an intuitive picture 
can straightforwardly be developed. Therefore, especially for young researchers, soft matter 
physics provides a promising opportunity in a meanwhile well-recognized environment. 
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